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1.  Introduction 


The  objective  of  this  research  program  is  to  determine  how  2-D  and  3-D 
disturbances  affect  the  mixing  of  species  in  two-stream  shear  layers  and  wakes.  The 
research  is  based  on  measurements  of  the  mixing  field  in  a  gravity-driven  liquid  mixing 
layer  facility.  A  major  part  of  this  effort  involves  the  quantification  of  the  mixing  field 
by  whole-field  measurements  using  laser  induced  fluorescence  (LIF)  diagnostics.  Mixing 
is  to  be  quantified  in  terms  of  the  concentration  field  of  a  passive  scalar,  using  its 
probability  density  function  (pdf)  and  the  derived  quantities  based  on  the  pdf.  While  the 
main  goal  here  is  to  investigate  basic  fluid  mechanical  phenomena,  it  is  expected  that  the 
knowledge  gained  will  be  useful  in  devising  active  and  passive  means  to  modify/control 
the  mixing  field  in  turbulent  shear  layers. 

This  final  report  provides  a  brief  summary  of  the  highlights  of  our  results  under  this 
grant.  Most  of  the  results  have  been  described  in  greater  detail  in  previous  annual  reports 
[1-2]  and  conference  presentations/papers  [3-4].  An  extensive  study  of  the  effects  of  2-D 
forcing  of  shear  layers  by  an  oscillating  airfoil  was  the  subject  of  the  recent  M.S.  thesis 
of  Mr.  Gregory  Katch  [5];  this  work  is  attached  in  the  Appendix. 

It  should  be  noted  that  the  work  carried  out  under  this  grant  is  being  continued  with 
the  support  of  AFOSR  Grant  No.  F49620-93- 1-0282,  which  itself  benefits  from  our  novel 
diagnostics  development  under  the  companion  AASERT  grant  (Grant  No.  F49620-92-J- 
0338).  Accomplishments  in  the  latter  two  research  grants  are  described  separately  in  their 
respective  annual  progress  reports. 


2.  Summary  of  2-D  Forcing  Results 

Analysis  and  comparison  of  data  from  two  sets  of  experiments  corresponding  to  two 
different  2-D  forcing  schemes  have  lead  to  a  unified  picture  of  the  effect  2-D  forcing  on 
the  shear  layer  mixing  field.  In  the  original  study  [6],  the  shear  layer  was  forced  by 
oscillating  the  high-speed  freestream  speed.  In  the  more  recent  experiments  [3,5],  2-D 
perturbations  were  introduced  into  the  shear  layer  by  a  2-D  oscillating  airfoil  placed 
downstream  of  the  splitter  plate  tip.  The  angle  of  attack  amplitude  A  and  the  frequency 
/  were  independently  controlled.  These  experiments  were  more  extensive  in  that  a  range  of 
forcing  amplitudes  was  considered  and  data  were  acquired  at  various  downstream  stations. 

The  growth  rate  characteristics  of  forced  shear  layers  have  been  described  in  the 
past  in  terms  of  a  non-dimensional  downstream  distance  x*  -  xXflU^,  where 
X  =  (f/j  - 1/2)  /  (^1  +  and  =  (I/j  + 1/2)/2.  The  flow  exhibits  an  enhanced  growth 
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rate  for  0  <  /  <  1,  the  region  1  <  /  <  2  refers  to  the  "frequency-locked"  region 
corresponding  to  the  passage  of  a  non-interacting  array  of  vortices  and  a  nearly  constant 
shear  layer  thickness,  and  jc*  >  2  describes  the  gradual  relaxation  to  the  growth 
characteristics  of  the  unforced  layer  through  various  vortex  interactions. 

Our  results  indicate  that  much  of  the  global  behavior  of  the  mixing  field  can  also 
be  described  in  terms  of  x*.  The  details  are  included  in  the  M.S.  thesis  attached  in  the 
Appendix.  To  highlight  the  role  of  jc*,  the  behavior  of  two  particular  aspects  of  mixing  - 
the  amount  of  mixed  fluid  in  the  center  of  the  shear  layer  and  the  total  amount  of  mixed 
fluid  across  the  layer  -  are  discussed  below. 

The  maximum  amount  of  mixed  fluid  is  found  in  the  middle  of  the  shear  layer  and 
can  be  characterized  in  terms  of  the  maximum  value  of  total  mixed-fluid  probability 
Pm(max)-  Figure  1  shows  that  the  amount  of  mixed  fluid  in  the  middle  of  a  forced  layer 
decreases  and  reaches  a  minimum  near  jc*  =  1.  The  total  amount  of  mixed  fluid  integrated 
across  the  layer,  quantified  in  terms  of  the  mixed-fluid  thickness  5^,  shows  a  significant 
increase.  Owing  to  a  comparable  increase  in  the  shear  layer  width  5j,  however,  the 
fraction  of  the  layer  width  occupied  with  mixed  fluid  (i.e.  5^/5i)  changes  only  modestly. 
These  results  are  illustrated  in  Figure  2. 

The  main  message  from  Figures  1  and  2  is  that  2-D  forcing  of  a  shear  layer  does 
not  necessarily  lead  to  enhanced  mixing.  In  fact,  the  region  of  enhanced  mixing  is  defined 
in  terms  of  x*.  In  the  range  jc*  <  1.5,  the  forced  layer  has  less  mixed  fluid  per  unit  width 
than  the  unforced  case.  The  maximum  amount  of  mixed  fluid  in  the  center  of  the  layer, 
Fm(max)’  the  total  amount  of  mixed  fluid  per  unit  width  reach  a  minimum  at 

about  JC*  =  1.  Figure  2  suggests  that  enhancement  of  5jn/5j  occurs  for  x*  >  1.5  which  can 
be  achieved  by  either  moving  downstream  (for  a  given  forcing  frequency),  or  by  increasing 
the  forcing  frequency  (for  a  fixed  downstream  station). 


3.  Summary  of  3-D  Forcing  Results 

Mixing  in  turbulent  shear  layers  is  closely  tied  to  the  generation  of  streamwise 
vorticity  as  part  of  the  process  leading  to  the  much  increased  interfacial  area  as  the  flow 
becomes  three-dimensional  and  turbulent.  The  investigation  into  3-D  forcing  is  motivated 
in  part  by  our  previous  results  that  2-D  forcing  alone  may  not  be  very  effective  for  mixing 
enhancement  and  control.  Our  research  to  date  has  concentrated  on  two  types  of  3-D 
effects;  those  generated  by  an  isolated  spanwise  disturbance  element  on  the  splitter  plate 
and  those  that  are  induced  by  the  side  walls  of  the  flow  facility. 
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3.1  3-D  Structure  and  Mixing  of  a  Forced  Wake  inside  a  Confined  Channel 

The  time  evolution  of  the  spanwise  flow  stmcture  and  the  mixing  field  in  a  low 
Reynolds  number  wake  (initial  Reg  =  100)  were  investigated.  New  information  was 
obtained  based  on  processing  a  large  data  set  corresponding  to  the  concentration  field 
^(x,y,z,f)  at  various  downstream  stations,  x,  over  the  cross-stream  plane  (i.e.  in  the  y-z 
plane)  as  a  function  of  time.  The  findings,  described  in  more  details  previously  [1,2,4], 
are  summarized  below. 

High  amplitude  2-D  forcing  of  the  low  Reynolds  number  wake  can  lead  to  a 
dramatic  increase  in  mixing.  The  streamwise  vorticity  in  this  flow  is  generated  by  the 
reorientation  and  stretching  of  the  main  2-D  spanwise  vorticity  as  a  result  of  the  shear 
region  in  the  side-wall  boundary  layer.  The  downstream  evolution  of  the  streamwise 
vorticity  produces  a  region  of  enhanced  mixing  which  is  confined  to  the  central  portion  of 
the  test  section  span,  with  the  regions  near  the  side  walls  practically  devoid  of  any  mixed 
fluid.  It  is  estimated  that  about  43%  of  the  test  section  cross  sectional  area  is  occupied 
by  mixed  fluid  at  the  farthest  downstream  measurement  station.  This  represents  an  order 
of  magnitude  increase  in  mixing  compared  with  the  unforced  case  at  the  same  downstream 
location.  It  can  be  envisioned  that  this  type  of  mixing  enhancement  may  play  a  role  in  an 
overall  strategy  for  mixing  control  in  confined  mixing  chambers.  Of  course,  we  still  have 
to  determine  the  range  of  flow  Reynolds  numbers  that  can  benefit  from  the  mixing 
enhancement  mechanism  identified  here. 

In  trying  to  understand  the  complex  3-D  structure  of  this  flow,  we  have  recently 
taken  advantage  of  the  flow  periodicity  in  constructing  the  detailed  "instantaneous"  3-D 
mixing  field  based  on  phase-averaged  LIF  images  of  the  flow.  A  3-D  reconstruction  of 
the  flow  stmcture  using  scientific  visualization  tools  are  being  used  to  aid  in  the 
understanding  of  the  flow  topology  near  the  wall.  These  developments,  which  started  under 
this  grant,  were  mostly  carried  out  with  support  from  grant  F49620-93- 1-0282.  Therefore, 
this  aspect  of  our  progress  will  be  more  appropriately  described  in  the  annual  report  of  the 
latter  grant. 


3.2  3-D  Forcing  of  a  Two-Stream  Shear  Laver  at  the  Splitter  Plate 

Among  the  various  types  of  splitter  plate  3-D  perturbation  mechanisms  investigated 
in  the  past,  we  decided  to  first  investigate  the  effects  on  the  mixing  field  generated  by  a 
single  disturbance  element.  A  cyhndrical  "peg"  was  placed  on  the  splitter  plate  similar  to 
that  used  by  Lasheras,  et  al.  [7]  in  a  low  Reynolds  number  shear  layer.  Mehta  &  Bell  [8] 
had  also  used  a  similar  device,  but  in  a  spanwise  periodic  arrangement.  The  peg  had  a 
diameter  D  =  1  mm,  a  height  H  =  3  mm,  and  was  placed  approximately  15  mm  upstream 
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of  the  splitter  plate  tip  at  the  mid-span  location  on  the  low-speed  side  of  the  splitter  plate. 
Except  for  the.  addition  of  the  peg,  the  flow  conditions  were  the  same  as  our  previous  2-D 
forced  shear  layers  (Uj  «  40  cm/s  and  U2  =  20  cm/s).  At  these  conditions,  the  height  of 
the  peg  exceeded  the  momentum  thickness  of  the  boundary  layer  in  which  it  was  placed. 
In  these  experiments,  the  shear  layer  is  forced  simultaneously  by  both  the  streamwise 
vorticity  due  to  the  peg  and  2-D  perturbations  by  oscillating  the  high-speed  freestream. 
Some  details  of  the  flow  structure  and  mixing  field  have  been  previously  reported  [1,2]. 
Three  particular  findings  constitute  our  main  results  and  they  are  summarized  below. 

The  addition  of  3-D  forcing  causes  a  spanwise  variation  in  the  total  amount  of 
mixed  fluid  characterized  by  a  mixed-fluid  thickness  variation  5^(z)  which  has  a  higher 
than  average  value  in  the  region  near  the  mid-span  location  (i.e.  the  location  of  the  3-D 
perturbation  source).  The  local  shear  layer  thickness  also  varies  in  the  spanwise  direction 
and,  in  this  case  thickness,  is  defined  by  8j(z).  The  spanwise  variation  of  the  amount  of 
mixed  fluid  per  unit  width  of  the  layer,  8^(z)/5i(z),  is  shown  in  Figure  3  for  various  forcing 
frequencies.  Note  that  the  effect  of  the  3-D  disturbance  on  the  splitter  plate  does  not 
remain  localized;  the  entire  span  is  affected.  The  well-defined  spanwise  wavy  pattern  is 
the  signature  of  the  counter-rotating  streamwise  vorticity  introduced  by  the  addition  of  3-D 
forcing.  As  the  forcing  frequency  increases,  the  mixed-fluid  fraction,  8^/81,  first  decreases 
and  then  increases  across  the  span.  This  trend  is  better  seen  when  the  span-averaged 
mixed-fluid  fraction  is  plotted  versus  the  forcing  frequency  /  and  the  corresponding  non- 
dimensional  distance  x  (see  Figure  4).  We  note  that  the  span-averaged  mixed-fluid 
fraction  in  this  flow  exhibits  the  same  dependence  on  x*  as  that  in  shear  layers  forced  by 
purely  2-D  disturbances.  The  highest  value  of  mixed-fluid  fraction  (30%  over  unforced) 
occurs  at  the  highest  forcing  frequency  of  /  =  16  Hz  which  corresponds  to  the  largest  x 
we  have  studied  to  date. 

The  presence  of  strong  three-dimensionality  in  the  case  of/=  16  Hz  leads  to  unique 
new  phenomena  which  are  revealed  when  the  spanwise  variation  of  the  pdf  of  the 
concentration  field  is  considered.  The  contour  plots  of  the  transverse  (y)  distribution  of  the 
pdf’s  at  three  spanwise  (z)  locations  are  shown  in  Figure  5.  The  most  striking  feature  is 
that  the  pdf’s  near  the  center-span  (z  =  0),  where  the  3-D  disturbance  source  is  located,  are 
non-uniform  across  the  layer  (i.e.  have  a  "marching"  behavior).  By  contrast,  the  pdf’s  at 
off-center  span  locations  (z  =  ±  12  mm)  exhibit  uniformity,  i.e.  a  "non-marching"  behavior. 
These  different  pdf  sttuctures  are  apparently  determined  by  the  relative  significance  of  2-D 
and  3-D  structures  in  the  transport  and  mixing  of  species,  a  finding  also  supported  by 
recent  computations  [9,10]. 

The  overall  (span-averaged)  pdf  of  concentration  field  for  the  highest  forcing 
frequency  displays  an  interesting  new  feature  for  this  flow;  see  Figure  6.  At  lower 
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frequencies,  the  predominant  mixed-fluid  concentration  does  not  vary  much  as  a  result  of 
forcing  in  agreement  with  previous  findings  [3,5,6].  At  /  =  16  Hz,  however,  there  is  a 
significant  shift  in  the  predominant  mixed-fluid  composition  to  lower  values  of 
concentration.  This  is  the  largest  modification  of  the  mixed  fluid  composition  we  have 
obtained  to  date.  We  do  not  currently  have  data  at  this  forcing  frequency  for  a  shear  layer 
perturbed  by  purely  2-D  disturbances.  Therefore,  we  cannot  yet  tell  whether  the  observed 
composition  shift  is  a  purely  2-D  effect  or  a  combined  2-D/3-D  effect.  Work  in  progress 
to  resolve  this  issue. 


4.  Work  in  Progress 

The  efforts  in  3-D  forcing  constitute  the  main  aspect  of  our  continuing  work.  We 
are  applying  scientific  visualization  methods  to  our  volumetric  data  sets  in  order  to 
understand  the  influence  of  forcing  on  the  flow  structure  and  mixing  field.  Forcing  affects 
both  the  spanwise  and  streamwise  structure  of  the  flow.  It  is  not  surprising,  therefore,  that 
we  are  finding  it  is  essential  that  the  flow  structure  be  viewed  in  3D  space  to  reveal  the 
intricate  interplay  between  streamwise  and  spanwise  structures. 

The  results  discussed  above  were  obtained  in  our  small  shear  layer  facility  where 
the  flow  Reynolds  number.  Re  =  AU6j/v,  is  limited  to  about  8,000.  Many  of  the  new 
findings  we  have  reported  need  to  be  investigated  at  higher  Reynolds  numbers.  A  new 
larger  shear  layer  facility  capable  of  Reynolds  numbers  at  least  an  order  of  magnitude 
higher  than  the  present  was  designed.  The  original  plans  for  performing  experiments  in 
this  facility  had  to  be  postponed  due  to  major  delays  in  the  delivery  of  the  facility.  The 
new  facility  has  finally  been  delivered  and  we  are  in  the  process  of  bringing  it  on-line. 
The  higher  Reynolds  number  experiments  will  be  carried  out  as  part  of  the  continuation 
grant  F49620-93- 1-0282. 
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Figure  1.  Amount  of  mixed  fluid  in  the  center  of  the  forced  shear  layer 
versus  x*;  forcing  amplitude  A  =  4  deg. 
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ABSTRACT 


MIXING  OF  SPECIES  IN  A  TWO-STREAM  SHEAR  LAYER 
FORCED  BY  AN  OSCILLATING  AIRFOIL 
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The  modifications  to  the  structure  and  the  mixing  field  in  a  shear  layer  forced  by 
an  oscillating  airfoil  were  documented  over  a  range  of  frequencies,  amplitudes  and 
streamwise  locations  using  laser  induced  fluorescence  in  the  passive  scalar  mode.  The 
shear  layer  structure  and  growth  were  observed  to  conform  to  the  well-established 
behavior  of  forced  shear  layers.  Results  on  the  mixing  field  show  a  variation  in  the 
amount  of  mixed  fluid  in  the  central  portion  of  the  forced  shear  layer  with  increasing 
downstream  distance  and  forcing  frequency,  an  increase  in  the  width  of  the  region 
occupied  by  mixed  fluid  and  an  increase  in  the  total  amount  of  mixed  fluid  in  the  layer. 
Only  for  the  highest  forcing  frequency  examined  in  this  study  was  there  an  increase  in 
the  fraction  of  the  layer  occupied  by  mixed  fluid  when  compared  to  the  unforced  layer. 
These  results  are  consistent  with  those  from  a  shear  layer  forced  by  a  different  method 
(oscillating  one  freestream).  Attempts  were  made  to  describe  the  downstream  evolution 
of  the  forced  mixing  layer  in  terms  of  a  universal  non-dimensional  curve.  The 
predominate  mixed  fluid  concentration  in  the  forced  shear  layer  showed  similar  behavior 
to  the  layer  forced  by  oscillating  one  freestream. 
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Chapter  1 

INTRODUCTION 


The  shear  layer  has  been  the  subject  of  a  large  number  of  studies  over  the  past 
few  decades.  One  reason  for  the  great  interest  in  this  flow  is  the  common  occurrence  of 
shear  layers  in  both  naturally  occurring  phenomena  and  in  industrial  processes.  Some 
familiar  examples  of  shear  layers  are  the  wind  blowing  over  a  body  of  water,  the  injection 
of  fuel  into  a  combustor  and  the  separated  flow  over  an  airfoil.  A  shear  layer  is 
generated  by  allowing  two  parallel  streams  of  fluid  moving  at  different  speeds  (initially 
separated  by  a  thin  interface)  to  come  together.  The  flow  formed  as  a  result  of  the  two 
fluids  merging  is  the  shear  layer.  A  schematic  of  the  geometry  of  a  plane  shear  layer  is 
shown  in  Figure  1. 

Consider  modelling  the  combustor  as  a  plane  shear  layer,  where  one  fluid  is  the 
fuel  and  the  other  fluid  is  the  oxidizer.  If,  for  example,  mixing  of  the  fuel  and  oxidizer 
in  the  combustor  could  be  optimized,  less  unburned  fuel  would  be  exhausted  from  the 
combustor,  resulting  in  a  more  efficient  combustion  process.  Like  the  combustor,  many 
other  industrial  applications  would  benefit  from  a  more  detailed  understanding  of  shear 
layer  dynamics  and  the  mixing  process  in  turbulent  shear  layers. 

Previous  studies  on  shear  layers  have  examined  many  different  aspects  of  natural 
and  forced  shear  layers.  In  the  following  sections  the  behavior  of  the  natural  layer,  the 
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forced  layer  and  mixing  in  the  shear  layer  will  be  discussed  as  related  to  this  study. 


1.1  The  natural  shear  layer 

It  was  first  noted  by  Brown  and  Roshko  (1971,  1974)  that  in  addition  to  the 
numerous  small  scales,  large  scale  structures  were  a  prominent  feature  of  the  natural  shear 
layer  over  a  wide  range  of  Reynolds  numbers.  Dimotakis  and  Brown  (1976)  showed  that 
the  large  scale  structures  seen  by  Brown  and  Roshko  persist  at  Reynolds  numbers  up  to 
3  X  10®,  based  on  high  speed  velocity  and  downstream  distance.  It  was  also  noted  that  the 
shear  layer  dynamics  were  governed  by  a  global  feedback  mechanism.  This  feedback 
mechanism  suggests  that  the  initiation  of  the  layer  at  the  splitter  plate  tip  is  coupled  to 
the  large  scale  structures  located  far  downstream. 

Growth  of  the  natural  shear  layer  was  observed  to  occur  as  the  result  of  the 
amalgamation  of  neighboring  large  scale  structures  by  Winant  and  Browand  (1974).  This 
process,  resulting  in  a  single  larger  structure,  termed  pairing,  results  (on  average)  in  the 
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linear  growth  of  the  natural  layer. 

Many  experimentalists  have  observed  large  scale,  two-dimensional  structures  in 
the  natural  shear  layer,  e.g.  Browand  and  Weidman  (1976)  and  Wygnanski,  Oster,  Fiedler 
and  Dziomba  (1979).  It  was  Miksad  (1972),  however,  who  first  observed  a  weaker 
longitudinal  (streamwise)  vortex  structure  in  the  (low  Reynolds  number)  natural  shear 
layer.  Later,  Konrad  (1976)  also  saw  these  streamwise  structures  at  much  higher 
Reynolds  numbers,  and  concluded  that  the  transition  to  turbulence  occurred  as  a  result  of 
the  formation  of  this  secondary  flow  structure.  Breidenthal  (1978,  1981)  also  observed 
the  streamwise  structures  and  suggested  that  this  structure  was,  in  fact,  pairs  of  counter 
rotating  vortices. 

1.2  The  forced  shear  layer 

Forcing  refers  to  the  application  of  a  periodic  disturbance  to  a  shear  layer  in  an 
attempt  to  modify  the  behavior  of  the  layer  by  exciting  the  natural  instabilities  of  the 
flow.  Much  work  has  been  done  on  understanding  forced  shear  layers  and  many  different 
methods  of  forcing  have  been  used.  Some  of  the  various  forcing  methods  previously  used 
include:  oscillating  one  or  both  of  the  freestream  speeds.  Ho  and  Huang  (1982),  Roberts 
and  Roshko  (1985)  and  Koochesfahani  and  M^'^Kinnon  (1991),  oscillating  the  splitter  plate 
tip,  Oster  and  Wygnanski  (1982),  forcing  by  acoustical  methods,  Fiedler  et  a/.  (1981)  and 
Zaman  and  Hussain  (1981),  spanwise  heating  on  the  surface  of  the  splitter  plate,  Nyggard 
and  Glezer  (1991),  and  by  oscillating  an  airfoil  downstream  of  the  splitter  plate, 
Koochesfahani  and  Dimotakis  (1989). 

Perhaps  the  most  important  observation  of  forced  shear  layers  is  the  control  over 
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the  growth  rate  and  the  structure  of  the  shear  layer  by  external  forcing.  It  was  reported 
by  Oster  and  Wygnanski  (1982)  and  Ho  and  Huang  (1982)  among  others,  that  low 
amplitude  two-dimensional  periodic  forcing  can  be  used,  to  some  extent,  to  control  the 
growth  and  structure  of  the  shear  layer.  The  effect  of  forcing  on  the  growth  of  the  shear 
layer  is  demonstrated  in  Figure  2,  reproduced  from  Browand  and  Ho  (1983).  In  this 
Figure  shear  layer  growth  is  quantified  by  the  local  momentum  thickness  of  the  shear 


Figure  2.  Effect  of  forcing  on  the  shear  layer  growth,  reprinted  from  Browand  and 
Ho  (1983). 
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layer  as  defined  by 
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In  this  figure  the  independent  axis  is  the  Wygnanski-Oster  parameter,  x*  =xXflU^, 
where  jc  is  the  distance  from  the  splitter  plate  tip,  X  =  {U ^ I {U ^+U ^  is  the  shear  layer 
growth  rate  parameter,  /  is  the  forcing  frequency  and  =  {U^+U^H.  This  parameter 
is  used  to  define  three  distinct  regions  of  the  evolution  of  the  forced  shear  layer. 

The  first  region,  called  the  enhanced  growth  region,  x*<  1,  exhibits  enhanced 
growth  rates  up  to  twice  that  of  the  natural  shear  layer.  The  next  region,  the  frequency 
locked  region,  1  <  x*  <2,  shows  reduced  growth  rates  when  compared  to  the  natural 
layer.  In  this  region  the  structures  are  observed  to  be  equally  spaced,  and  have  passage 
frequency  equal  to  the  forcing  frequency.  In  the  third  region,  x*>  2,  the  growth  rate 
relaxes  to  the  natural  growth  rate. 

1.3  Mixing  in  the  shear  layer 

A  majority  of  the  experimental  and  computational  studies  to  date  have  focused  on 
the  effects  of  forcing  on  the  momentum  transport  properties  of  the  shear  layer  such  as  the 
layer  growth  rate,  the  resulting  velocity  and  vorticity  fields  and  their  fluctuating 
components  and  the  Reynolds  stresses.  By  contrast,  there  have  been  few  detailed  studies 
devoted  specifically  to  the  behavior  of  the  scalar  mixing  field  in  forced  shear  layers. 
Throughout  this  work,  mixing  will  refer  to  the  mixing  of  species,  the  kind  of  mixing  that 
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is  necessary  for  chemical  reaction  and  combustion. 

Konrad  (1976)  and  Breidenthal  (1978,  1981)  were  among  the  first  to  consider 
molecular  mixing  in  the  shear  layer.  Konrad  (1976)  found  in  a  non-reacting  gaseous 
layer,  a  rapid  increase  in  the  amount  of  mixed  fluid  at  some  distance  downstream  of  the 
splitter  plate.  The  increase  in  mixing  was  attributed  to  increased  interfacial  area  between 
the  two  fluids  resulting  from  the  development  of  small-scale  motion  which  were  generated 
during  transition  to  turbulence.  Breidenthal  (1978,  1981)  and  Koochesfahani  and 
Dimotakis  (1986)  also  found  this  dramatic  increase  in  mixing,  in  chemically  reacting 
liquid  layers,  and  referred  to  this  as  the  mixing  transition.  The  mixing  transition  is  shown 
in  Figure  3,  reproduced  from  Roshko  (1990).  The  velocity  ratio  r  =  C/j  IU2  in  Figure  3 
is  2.6. 
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Figure  3.  The  mixing  transition,  reprinted  from  Roshko  (1990). 
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The  velocity  ratio,  r  =  t/i  IU2 ,  has  been  shown  to  have  an  effect  on  the  mixing 
transition,  see  Figure  4,  reproduced  from  Breidenthal  (1978).  In  Figure  4  the  velocity 
ratio  is  defined  as  r  =  U^IU^ ,  the  inverse  of  the  notation  used  in  this  work.  In  the  current 
work  the  velocity  ratio  is  r  =  t/j  IU2  =  2. 


Figure  4.  Effect  of  the  velocity  ration  on  the  mixing  transition,  reprinted 
from  Breidenthal  (1978). 


Roberts  (1984,  1985)  observed  the  amount  of  chemical  product  in  a  liquid  shear 
layer  using  chemically  reacting  techniques.  It  was  found  that  the  effect  of  two- 
dimensional  forcing  can  significantly  alter  the  amount  of  chemical  product  in  the  layer. 
In  pre-transitional  flows  there  were  large  increases  in  the  amount  of  chemical  product  in 
the  frequency  locked  region.  In  post-transitional  flows  increases  were  observed  only  in 
the  very  early  stages  of  the  enhanced  growth  region.  In  the  frequency  locked  region  and 
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beyond  there  was  actually  less  chemical  product  in  the  forced  layer  than  in  the  natural 
layer.  The  technique  of  Laser  Induced  Fluorescence  (LIF)  introduced  by  Koochesfahani 
and  Dimotakis  (1986),  was  used  to  measure  the  chemical  product  in  reacting  flows  and 
the  composition  of  mixed  fluid  in  non-reacting  flows. 

The  effects  of  periodic  oscillation  of  the  high-speed  freestream  on  the  composition 
of  mixed  fluid  in  a  non-reacting  shear  layer  were  reported  by  Koochesfahani  and 
M“‘^Kinnon  (1991).  It  was  found  that  while  the  total  amount  of  mixed  fluid  integrated 
across  the  layer  had  increased,  the  amount  of  mixed  fluid  per  unit  width  of  the  layer  had, 
in  fact,  remained  nearly  constant.  This  suggests  that  the  larger  growth  rate  of  a  shear 
layer  forced  by  a  two-dimensional  disturbance  does  not  necessarily  lead  to  a  more 
efficient  mixer.  In  this  work  efficiency,  as  related  to  the  mixing  of  species  in  a  shear 
layer,  is  defined  as  the  amount  of  mixed  fluid  per  unit  width  of  the  layer.  It  was  also 
noted  that  forcing  shifted  the  predominant  mixed  fluid  concentration  to  higher  values,  i.e. 
larger  proportion  of  high-speed  to  low-speed  fluid. 

1.4  Objectives 

As  stated  earlier,  there  are  many  different  ways  to  perturb  a  shear  layer.  All  of 
the  methods  previously  discussed  impose  disturbances  on  the  shear  layer  by  different 
mechanisms,  each  entering  the  flow  in  a  different  way.  Surprisingly,  the  net  result  on  the 
growth  rate  and  structure  of  the  shear  layer  in  response  to  forcing  is  very  similar  for  all 
methods.  The  goal  of  the  present  work  is  to  investigate  the  influence  of  the  forcing 
mechanism  on  the  mixing  of  species  in  a  two  stream  shear  layer. 

In  an  earlier  investigation,  the  effects  of  the  periodic  oscillation  of  the  high-speed 
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freestream  on  the  composition  of  mixed  fluid  in  a  plane  shear  layer  were  reported  by 
Koochesfahani  and  M‘“'Kinnon  (1991).  The  results  from  the  oscillating  freestream 
experiments  are  compared  to  a  shear  layer  forced  by  oscillating  an  airfoil  located 
downstream  of  the  splitter  plate  tip.  In  addition  to  comparing  the  methods  of  forcing,  the 
effect  of  the  forcing  amplitude  and  streamwise  evolution  on  the  total  amount  of  mixed 
fluid  in  the  layer,  the  amount  of  mixed  fluid  in  the  central  portion  of  the  layer  and  the 
composition  of  mixed  fluid  will  also  be  examined. 


Chapter  2 

EXPERIMENTAL  FACILITY  AND  INSTRUMENTATION 


All  of  the  data  presented  in  this  work  were  acquired  over  a  two  day  period,  in  two 
sets  of  experiments.  Each  of  the  two  experiments  consisted  of  ten  forcing  conditions  and 
covered  a  different  streamwise  view  of  the  test  section.  In  this  chapter  the  shear  layer 
facility,  forcing  mechanism  and  diagnostics  are  described. 

2.1  Experimental  facility 

The  experiments  were  performed  in  a  gravity-driven  water  shear  layer  apparatus, 
a  schematic  of  the  facility  is  shown  in  Figure  5.  For  reference,  the  inside  dimensions  of 
the  test  section  were  4  (height)  x  8  (width)  x  36  (length)  cm.  The  system  was  fed  from 
two  210  liter  reservoirs,  one  for  each  the  high-speed  and  low-speed  sides  of  the  shear 
layer.  Water  was  pumped  from  the  reservoirs  to  constant  head  tanks,  located 
approximately  eight  feet  above  the  test  section.  Flow  through  the  test  section  was 
regulated  by  three  valves;  one  valve  in  each  of  the  supply  lines  upstream  of  the 
contraction,  and  a  third  valve  downstream  of  the  test  section. 

For  these  experiments  the  freestream  velocities  were  set  to  =  40  cm/s  and  U2 
=  20  cm/s,  yielding  a  velocity  ratio  of,  r  =  2.  In  this  work  the  flow  was  examined  over 
the  range  12.5  <  x  <  20.0  cm,  corresponding  to  0.28  <  x  *  <  1.78.  The  Reynolds  number. 
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Constant  Head  Tanks 


Re.  =  AU  SJv  ,  based  on  velocity  difference  and  shear  layer  width,  ranges  from  4400 
<  Re.  <  6000.  This  range  of  Reynolds  numbers  corresponds  to  the  later  stages  of  the 
mixing  transition,  see  Breidenthal  (1981)  and  Figure  4.  It  is  not  clear  however,  if  the 
presence  of  the  airfoil  affects  the  mixing  transition.  Due  to  the  constraints  of  the  flow 
facility,  this  is  nearly  the  highest  Reynolds  number  currently  attainable.  The  natural  shear 
layer  roll  up  frequency  f„ ,  at  the  splitter  plate  tip  was  approximately  27  Hz. 

An  important  consideration  to  flow  quality  was  the  development  of  the  boundary 
layer  on  both  the  upper  and  lower  surfaces  of  the  splitter  plate.  A  series  of  screens  and 
straws  were  used  to  create  a  uniform  flow  through  the  contraction.  Additionally,  extreme 
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care  was  taken  during  the  filling  process  to  expel  all  air  from  the  system.  Any  air  trapped 
under  the  splitter  plate  or  in  the  straws  may  affect  the  boundary  layer  growth  and  the 
resulting  shear  layer. 

2.2  Forcing  mechanism 

Perturbations  were  introduced  into  the  shear  layer  by  a  15%  thick,  2-D  airfoil  of 
chord  C  =  2  cm,  spanning  the  entire  width  of  the  test  section  b  =  S  cm,  oscillating 
sinusoidally  in  pitch  about  the  quarter-chord  with  angle  of  attack  amplitude  A  and 
frequency  /.  The  airfoil  was  rotated  about  its  quarter-chord  point,  located  midway 


k -  6.5  cm  - >\ 


Figure  6.  Schematic  of  the  forcing  arrangement. 
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between  the  upper  and  lower  surfaces  of  the  test  section  and  6.5  cm  downstream  of  the 
splitter  plate  tip.  A  schematic  of  the  forcing  mechanism  is  shown  in  Figure  6.  The 
placement  of  the  airfoil  was  determined  by  the  physical  constraints  of  the  airfoil 
assembly.  Movement  of  the  airfoil  was  achieved  by  an  electro-magnetic  coil  vibrator 
(Vibrations  Test  Systems  VG  100-8).  To  improve  the  sinusoidal  motions  of  the  airfoil 
a  P-D  feedback  control  system  was  used  to  match  the  airfoil  position  to  a  command  input. 

The  angle  of  attack  of  the  airfoil  varied  between  -A  and  -i-A,  with  a  mean  angle 
of  attack  approximately  zero.  In  this  work  an  unforced  case  and  nine  forced  cases  were 
examined.  In  the  unforced  case  the  airfoil  was  stationery  at  zero  degrees  angle  of  attack. 
The  forced  cases  included  three  frequencies  at  three  amplitudes  each;  the  combinations 
examined  are  shown  in  Table  1. 


Table  1.  Forcing  frequency-amplitude  combinations. 


/(Hz) 

A  (degrees) 

2 

2,  4,  8 

4 

2,  4,6 

8 

2,  4,6 

2.3  Diagnostics 

The  mixing  field  was  measured  using  Laser  Induced  Fluorescence  (LIF)  in  the 
passive  scalar  (non-reacting)  mode.  This  technique  utilizes  a  passive  scalar  containment, 
such  as  a  fluorescent  dye  that  is  initially  mixed  with  one  of  the  freestream  fluids.  In  this 
case  the  low-speed  fluid  was  marked  with  the  fluorescent  dye  disodium  fluorescein.  The 
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dye  is  then  diluted  as  a  result  of  mixing  with  fluid  from  the  high-speed  stream.  By 
recording  the  fluorescence  intensity  issuing  from  a  sample  volume  as  a  function  of  time, 
a  quantitative  measurement  of  the  dye  concentration,  and  therefore  the  relative 
concentration  of  high-speed  to  low-speed  fluid  within  the  volume  may  be  made.  The 
instantaneous  dye  concentration  Q  in  the  sample  volume  is  defined  by 


where  C ,  is  the  freestream  (low-speed)  dye  concentration,  Vj  and  V2  are  the  volume  of 
fluid  from  the  high-  and  low-speed  freestreams,  respectively,  within  the  sample  volume. 
Throughout  this  work  the  concentration  ^  will  be  normalized  as 


^  = 


In  terms  of  the  high-speed  volume  fraction,  concentration  may  be  written  as 


('^1  +  ^2)' 


Therefore,  pure  fluid  from  the  low-  and  high-  speed  sides  have  concentrations  ^  =  0  and 
^  =  1,  respectively. 

The  LIF  measurements  were  carried  out  over  a  plane  defined  by  a  laser  sheet, 
approximately  0.5  mm  thick,  aligned  along  the  flow  direction  at  the  shear  layer  mid-span 
location.  The  laser  sheet  was  formed  by  passing  a  3  watt  beam  from  a  Lexel  model  95 
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(30  frames/sec)  with  an  exposure  time  of  2  msec.  Images  were  then  digitized  to  8  bits 
and  captured  onto  hard  disk  in  real  time  by  an  image  acquisition  system  (Recognition 
Concepts,  Inc.,  TRAPIX-5500).  For  each  of  the  ten  forcing  conditions,  256  consecutive 
LIF  images  (512  fields)  were  acquired,  corresponding  to  32,  64,  and  128  large  scale 
structures  passing  the  field  of  view  for  forcing  frequencies  of/ =  2,  4,  8  Hz,  respectively. 

The  4  cm  test  section  height  was  imaged  onto  200  of  480  vertical  pixels;  512 
pixels  correspond  to  approximately  10  cm  in  the  streamwise  direction.  The  resulting 
resolution  is  approximately  200  x  200  pm.  It  should  be  noted  that  at  this  resolution  the 
smallest  expected  diffusion  (Batchelor)  scales,  are  not  resolved.  According  to  Miller 
and  Dimotakis  (1991),  the  smallest  viscous  scale,  and  the  Kolmogorov  scale,  « 
S^Re  are  related  by  X^  ~  25  A^.  The  smallest  diffusion  scale  A^  is  another  factor 
of  Sc  smaller  than  the  smallest  viscous  scale,  resulting  in  A^  «  A^  for  water  (Sc  « 
600).  Over  the  range  of  parameters  examined  in  this  study  the  smallest  diffusion  scales, 
Xjj,  are  estimated  to  be  about  41  pm,  approximately  5  times  smaller  than  the  spatial 
resolution  of  the  current  memeasurement. 

A  curious  limitation  of  this  technique  is  that  the  measurements  always  over¬ 
estimate  the  amount  of  mixed  fluid.  Unmixed  low-  and  high-speed  fluids  may 
simultaneously  exist  within  the  measurement  resolution  and  appear  as  if  they  were 
actually  mixed  at  that  ratio.  The  results  presented  in  this  work  are,  therefore,  not  intended 
to  provide  an  absolute  measure  of  the  extent  of  molecular  mixing.  The  primary  focus  of 
these  results  is  to  characterize  the  relative  changes  in  the  mixing  field  in  a  forced  layer 
compared  to  the  unforced  layer.  Such  comparisons  are  believed  to  be  warranted  since  the 
resolution,  compared  to  the  smallest  diffusion  scale,  is  the  same  for  all  of  the 
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measurements  presented  in  this  work. 


2.4  Data  reduction 

The  first,  and  perhaps  the  most  time  consuming  step  in  the  data  reduction  process, 
is  the  normalization  of  the  images.  The  images  were  first  processed  to  remove  the 
Gaussian  distribution  of  the  laser  sheet  and  nonuniform  pixel  response  of  the  CCD 
camera.  A  calibration  was  performed  at  the  start  of  the  experiment  to  capture  these 
effects.  During  the  calibration  the  test  section  was  flooded  with  low-speed  fluid  mixed 
with  the  dye,  images  were  recorded  and  then  averaged  to  produce  a  calibration  image. 
The  nonuniformities  in  the  pre-processed  images  were  then  removed  on  a  pixel  by  pixel 
basis  according  to 


where  is  the  corrected  intensity,  I  is  the  instantaneous  image  intensity  and  is  the 
intensity  from  the  calibration  image.  Concentration  may  then  be  computed  as 


where  I  is  the  corrected  freestream  dye  intensity.  This  expression  is  simply  a 

COtTo 

conversion  from  the  measured  quantity,  fluorescence  intensity  to  dye  concentration.  The 
absorption  coefficient  a  =  e^C^,  where  Cq  is  the  dye  molar  absorption  coefficient  and 
is  the  molar  dye  concentration,  was  measured  by  Koochesfahani  (1984)  to  be  0.157  cm  * 
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at  a  dye  concentration  of  Q  =  10"^M.  Attenuation  of  the  fluorescence  intensity  at  a  given 
depth  y  is  then  given  by  /  =  IqC  ,  where  /q  is  the  initial  intensity.  In  the  current  work, 
the  dye  concentration  is  approximately  2.5xlO  ’M;  a  factor  of  40  smaller  than  in  the  work 
of  Koochesfahani  (1984).  In  the  event  that  the  test  section  is  completely  filled  with  dye, 
the  worst  possible  case,  the  fluorescence  intensity  at  the  bottom  of  the  test  section  would 
have  been  attenuated  by  only  1.6%.  The  effect  of  attenuation  is  much  smaller  than  the 
overall  signal-to-noise  ratio  of  the  experiment  e  =  4.5%. 

From  the  normalized  images  the  probability  density  function  (pdf)  of  the 
concentration  field  can  be  constructed  at  any  of  the  measurement  points  in  the  x-y  plane. 
The  pdf  of  concentration  ^  as  a  function  of  position  y  in  the  layer,  written  p{i,y),  was 
computed  at  six  different  streamwise  locations.  The  pdf  p  y)  is  the  basis  for  all  of 
the  computed  quantities  in  this  study.  The  ideal  pdf  would  have  delta  functions  at  ^  = 
0  and  1  representing  pure  fluid  from  the  low-  and  high-speed  freestreams,  respectively. 
In  practice,  the  delta  functions  have  a  finite  width  e  which  is  the  result  of  the  overall 
signal-to-noise  ratio  of  the  experiment.  In  this  work  e  =  0.045,  therefore  concentrations 
in  the  range  0  <  ^  <  e  and  1-e  <  ^  <  I  are  considered  to  be  pure  unmixed  fluid  from  the 
low-  and  high-speed  freestreams,  respectively.  Likewise,  concentrations  within  the  range 
e  <  ^  <  l-£  are  considered  mixed  fluid. 

The  first  quantity  computed  from  p  y)  is  the  average  concentration  ^  which 
is  found  by  integrating  over  all  concentrations  at  a  given  position,  as  defined  by 

1 

^(y)  =  f  ^piLy)d^. 

0 
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To  characterize  the  effect  of  forcing  on  the  mixing  field  the  probability  of  occurrence  of 
unmixed  fluid  from  either  the  low-  or  high-speed  freestreams  can  be  computed  by 

e 

Po(y)  =  / 

0 

1 

Pliy)  =  /  p(^,y)dt 

1-e 

Similarly,  the  probability  of  occurrence  of  mixed  fluid  is  computed  by  integrating  over 
the  range  of  concentrations  defining  mixed  fluid  as  given  by 

1-8 

Pm(y)  =  f  p(^’y)d^- 

e 

It  should  be  noted  that  p„  is  the  probability  of  occurrence  of  mixed  fluid  of  all 
concentrations;  a  value  of  less  than  one  implies  the  presence  of  unmixed  fluid. 

The  shear  layer  width  may  be  computed  based  on  the  total  mixed  fluid 
probability  curve  at  a  given  downstream  location.  The  local  width  is  defined  in  a  manner 
similar  to  boundary  layer  thickness.  Recall,  the  boundary  layer  thickness  is  the  height 
above  a  surface  where  the  local  velocity  is  99%  of  the  free  stream  velocity.  The  shear 
layer  width  Sj  is  similarly  defined  as  the  distance  between  the  two  points  where  the 
probability  of  mixed  fluid  is  1%  of  the  maximum  value  of  (y).  Koochesfahani  and 
Dimotakis  (1986)  reported  that  the  visual  width  of  the  shear  layer  is  approximately 
equal  to  bj. 

A  measure  of  the  total  amount  of  mixed  fluid  (at  all  concentrations)  in  the  layer. 
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the  mixed  fluid  thickness,  is  defined  by 

+h  1-  e  +h 

f  /  P{i,y)dldy=  f  p^(y)dy. 

-h  E  -h 

Note  that  the  mixed  fluid  thickness  is  simply  the  area  under  the  total  mixed  fluid 
probability  curve.  The  quantity  termed  the  mixed  fluid  fraction,  is  used  to 

quantify  the  amount  of  mixed  fluid  per  unit  width  of  the  shear  layer. 

It  has  been  shown  by  Koochesfahani  and  M"‘=Kinnon  (1991)  that  the  distribution 
of  the  mixed  fluid  composition  is  nearly  invariant  across  the  entire  width  of  the  layer,  and 
therefore  may  be  characterized  by  an  average  pdf.  The  average  pdf  is  defined  by 


where  h  is  the  test  section  half  width. 

A  common  measure  of  shear  layer  growth  is  the  momentum  thickness  defined  by 


e  = 


-00 


{U,-U){U-U,) 


In  this  work  no  velocity  measurements  were  made,  therefore  the  momentum  thickness 
cannot  be  calculated.  Another  integral  thickness,  called  the  concentration  thickness,  may 
be  defined  in  a  manner  analogous  to  the  momentum  thickness,  based  on  the  average 


concentration,  defined  as 
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+  00 


-00 


Recall  that  the  subscripts  1  and  2  in  the  above  equation  refer  to  the  average  concentration 
of  the  high-  and  low-speed  freestreams,  respectively.  The  average  concentration  of  the 
high-speed  freestream  =  1  and  for  the  low-speed  freestream  =  0.  The 
concentration  thickness  is  then  reduced  to 

+  00 

-00 

The  relationship  between  the  shear  layer  width  8^  and  the  concentration  thickness  0^  will 


be  examined  in  section  3.4. 


Chapter  3 

RESULTS  AND  DISCUSSION 


3.1  Flow  visualization  results 

Figures  8-12  illustrate  the  influence  of  the  forcing  frequency  and  amplitude  on 
the  structure  and  growth  of  the  shear  layer.  In  all  of  the  images  the  flow  is  from  right 
to  left  with  the  high-speed  stream  on  top.  The  time  between  images  At  is  1/30  sec,  (every 
other  field)  with  time  increasing  from  top  to  bottom.  Each  image,  composed  of  90,000 
pixels  (450  in  x  and  200  in  y),  spans  the  entire  height  of  the  test  section  and  the  range 
12  <  X  <  22  cm  (right  to  left).  After  processing,  the  images  were  pseudo-colored  to 
simulate  a  chemical  reaction.  Pure  fluid  from  either  of  the  freestreams  was  assigned  to 
black.  Fluid  composed  of  mainly  low-speed  fluid,  ^  ~  s,  were  assigned  shades  of  blue, 
likewise  nearly  pure  high-speed  fluid,  ^  =  1-e,  were  colored  shades  of  red  (recall  pure 
low-speed  fluid  has  concentration  0  <  ^  <  e,  high-speed  fluid  l-e  <  ^  <  1). 

Two  different  time  series  of  images  from  the  unforced  layer  are  presented  in 
Figure  8.  These  clearly  show  the  random  nature  of  the  unforced  layer.  In  Figure  8  (a) 
a  fairly  well  defined  large  structure  may  be  seen  passing  the  field  of  view,  while  some 
time  later  in  Figure  8  (b),  there  are  no  well  defined  large  structures  apparent.  It  should 
be  noted  that  the  unforced  shear  layer  described  here  is  not  the  same  as  the  natural  layer 
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due  to  the  presence  of  the  airfoil.  Differences  between  the  unforced  and  natural  shear 
layers  will  be  discussed  in  Section  3.6. 

Figures  9-12  display  the  effect  of  the  forcing  amplitude  A  on  the  structure  and 
growth  of  the  shear  layer  forced  at/=  2,  4  and  8  Hz,  respectively.  For  the  case  forced 
at  /  =  2  Hz,  A  =  2°,  see  Figure  9  (a),  there  is  very  little  if  any  visible  effect  on  the 
structure  and  growth  rate  of  the  layer,  when  compared  to  the  unforced  case.  With 
increasing  amplitude.  Figures  9  (b)  and  (c),  the  layer  growth  rate  increases  resulting  in 
larger  vortical  structures  and  shear  layer  width  at  a  given  downstream  location  x.  As  a 
result  of  the  larger  layer  width  mixed-fluid  can  also  be  seen  over  a  larger  portion  of  the 
test  section.  Quantitative  descriptions  of  the  shear  layer  width  and  the  amount  of  mixed 
fluid  will  be  given  in  section  3.3.  There  are  some  periodic  structures  in  the  flow, 
however  they  are  not  very  well  defined.  The  location  where  the  structures  are  fully 
formed  is  the  center  of  the  frequency  locked  region,  x*  =  1.5;  for  /  =  2  Hz  this 
corresponds  to  x  =  68  cm,  which  lies  well  beyond  the  end  of  the  test  section.  See  section 
1.2  for  description  of  the  nondimensional  streamwise  coordinate  x*.  Moving  farther 
downstream,  if  possible,  would  not  allow  the  structures  to  further  develop  because  of  the 
limiting  height  of  the  test  section.  The  lower  wall  of  the  test  section  may  be  seen  to 
influence  the  structures  in  Figure  9  (c).  In  the  images,  the  nondimensional  streamwise 
coordinate  spans  the  range,  0.26  <  x  *<  0.52,  which  lies  well  within  the  enhanced  growth 
region. 

In  Figure  10  the  structure  of  the  shear  layer  is  shown  as  a  function  of  the  forcing 
amplitude  A  when  forced  at/  =  4  Hz.  At  this  frequency  the  nondimensional  streamwise 
coordinate  spans  the  end  of  the  enhanced  growth  rate  region,  0.52  <x*<  1.02.  Again, 


I 
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Figure  9.  Effect  of  the  forcing  amplitude  on  the  structure  of  the  shear  layer  forced  at 
/  =  2  Hz,  flow  is  from  right  to  left,  At  =  1/30  sec. 


Figure  10.  Effect  of  the  forcing  amplitude  on  the  structure  of  the  shear  layer  forced 
at/=  4  Hz,  flow  is  from  right  to  left,  At  =  1/30  sec. 


Figure  11.  Effect  of  the  forcing  amplitude  on  the  structure  of  the  shear  layer  forced 
at/=  8  Hz,  flow  is  from  right  to  left,  At  =  1/30  sec. 
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it  may  be  seen  that  as  the  amplitude  increases  the  layer  growth  rate  and  the  amount  of 
mixed  fluid  within  the  layer  also  increases.  In  this  case  the  structures  are  more  well 
defined  than  the  case  forced  at  2  Hz.  At  the  largest  amplitude  the  structures  are,  again, 
approaching  the  width  of  the  test  section. 

The  cases  forced  at/  =  8  Hz,  shown  in  Figure  11,  span  the  frequency  locked 
region  (1.04  <x*<  2.04)  characterized  by  an  array  of  non-interacting,  equally  spaced 
vortices  and  reduced  growth  rate.  Even  at  the  smallest  amplitude  A  =  2°,  the  structures 
are  very  well  defined.  As  in  the  previous  cases,  the  effect  of  amplitude  is  shown  to  result 
in  an  increase  in  the  layer  width  and  the  amount  of  mixed  fluid  within  the  layer.  In  all 
of  these  cases,  the  increase  in  shear  layer  growth  rate,  in  response  to  low  frequency 
forcing  (relative  to  the  natural  frequency,  fg  =  27  Hz),  culminating  in  the  formation  of 
large  vortical  structures  is  consistent  with  previous  shear  layer  studies  eg.  Fiedler  et  al. 
(1981),  Oster  and  Wygnanski  (1982),  Ho  and  Huang  (1982),  Koochesfahani  and 
Dimotakis  (1989).  For  completeness.  Figure  12  is  included  to  show  the  effect  of  the 
forcing  frequency  on  the  structure  of  the  shear  layer  when  forced  at  ?l  fixed  amplitude  of 
A  =  4°.  Note  as  the  frequency  is  increased,  the  excursions  of  low-speed  fluid  to  the  upper 
side,  and  high-speed  fluid  to  the  lower  side  of  the  shear  layer  becomes  deeper.  This 
excursion  of  fluid  is  often  referred  to  as  an  entrainment  tongue. 

3.2  Transverse  profiles  of  Po  p^ 

The  quantities  computed  from  the  probability  density  function  of  the  high-speed 
fluid  volume  fraction are  presented  in  this  section.  The  plots  are  divided  into  two 
sections.  The  first  section  will  present  the  evolution  of  the  various  computed  quantities 
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as  a  function  of  streamwise  location  x,  at  a  fixed  forcing  frequency  and  amplitude.  The 
second  will  show  the  amplitude  dependence  at  a  fixed  downstream  location  x.  Data  were 
computed  at  6  uniformly  spaced  streamwise  locations  spanning  the  range  0.28  <  x*  <  1.78 
(12.5  <  X  <  20.0  cm). 

3.2.1  Streamwise  dependence 

Transverse  profiles  of  the  average  concentration  the  total  mixed  fluid  probability 
(y)  and  unmixed  fluid  probabilities,  Po(y)  and  p^  (y)  for  the  unforced  shear  layer  are 
shown  in  Figure  13.  In  general,  the  results  for  the  unforced  layer  are  very  similar  in 
shape  to  the  results  of  the  natural  layer  reported  by  Koochesfahani  and  Dimotakis  (1986) 
and  Koochesfahani  and  M^'Kinnon  (1991).  As  previously  discussed,  differences  in  the 
shear  layer  width  and  the  amount  of  mixed  fluid  between  the  unforced  and  natural  layers 
will  be  quantified  in  section  3.6. 

The  most  obvious  feature  of  the  unforced  layer  is  the  increase  in  the  layer  width 
with  increasing  downstream  distance  as  seen  in  the  total  mixed  fluid  probability  curves 
displayed  in  Figure  13  (b).  Increasing  shear  layer  width  may  also  be  inferred  from 
outward  movement  of  the  probability  of  pure  low-  and  high-speed  curves,  see  Figure  13 
(c)  and  (d).  A  feature  unique  to  the  unforced  layer  is  the  nearly  constant  value  of  the 
peak  of  the  mixed  fluid  probability  curve  0.95  with  increasing  downstream 

distance.  Recall  that  a  value  of  p„  less  than  unity  implies  the  presence  of  unmixed  fluid; 
therefore  roughly  5%  of  the  fluid  in  the  center  of  the  unforced  layer  is  unmixed.  It  will 
be  shown  later,  that  there  can  be  a  significant  streamwise  variation  in  in  the 
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(a)  Average  concentration 


(b)  Total  mixed  fluid  probability 


(c)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 

o  X  =  12.5,  □  X  =  14.0,  V  X  =  15.5,  •  x  =  17.0,  ■  x  =  18.5,  a  x  =  20.0  cm 


Figure  13.  Transverse  profiles  for  the  unforced  shear  layer  as  a  function  of  the 
downstream  distance  x. 
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(b)  Total  mixed  fluid  probability 
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(c)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 

o  jc*  =  0.28,  □  jc*  =  0.31,  V  jc*  =  0.34,  •  ;c*  =  0.37,  ■  x*  =  0.41,  a  x*  =  0.44 

Figure  14.  Transverse  profiles  for  the  shear  layer  forced  at  /=  2  Hz,  A  =  2°  as  a 
function  of  the  nondimensional  downstream  distance  x*. 
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(c)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 

o  jc*  =  0.56,  o  X*  =  0.62,  v  /  =  0.68,  •  x  =  0.75,  ■  x*  =  0.82,  a  ;c*  =  0.88 


Figure  17. 


Transverse  profiles  for  the  shear  layer  forced  at  /  =  4  Hz,  A  =  2°  as  a 
function  of  the  nondimensional  downstream  distance  x*. 


(a)  Average  concentration 


(b)  Total  mixed  fluid  probability 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


o  X*  =  0.56,  a  X*  =  0.62,  v  x* 


=  0.68,  •  X*  =  0.75,  mx*  =  0.82,  a  x*  =  0.88 


Figure  18.  Transverse  profiles  for  the  shear  layer  forced  at  /  =  4  Hz,  A  =  4°  as  a 
function  of  the  nondimensional  downstream  distance  x*. 
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(c)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 

o  X*  =  1.11,  □  X  =  1.24,  V  X  =  1.38,  •  X*  =  1.51,  ■  x*  =  1.64,  a  x*  =  1.78 

Figure  21.  Transverse  profiles  for  the  shear  layer  forced  at  /  =  8  Hz,  A  =  4°  as  a 
function  of  the  nondimensional  downstream  distance  x*. 
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(c)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 

o  X*  =  1.11,  a  X*  =  1.24,  V  X*  =  1.38,  •  x*  =  1.51,  ■  x*  =  1.64,  a  x*  =  1.78 

Figure  21.  Transverse  profiles  for  the  shear  layer  forced  at  /  =  8  Hz,  A  =  6°  as  a 
function  of  the  nondimensional  downstream  distance  x*. 
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Results  for  the  forced  cases  are  displayed  in  Figures  14  -  22.  A  noteworthy 
feature  of  the  forced  shear  layer  is  the  severe  distortion  of  the  average  concentration 
profiles  |  (y).  The  development  of  a  'flat  region'  in  the  center  of  the  ^  (y)  profile  is 
characteristic  of  the  passage  of  well  mixed  vortex  cores,  previously  noted  by  Wygnanski, 
Oster  and  Fiedler  (1979)  and  Koochesfahani  and  M^'^Kinnon  (1991).  This  'flat  region'  is 
seen  developing  in  the  cases  forced  at/  =  4  Hz  and  is  fully  developed  in  the  8  Hz  cases. 

Similar  to  the  unforced  layer,  the  forced  cases  all  show  increasing  shear  layer 
width  as  downstream  distance  increases.  Correspondingly,  mixed  fluid  is  found  over  a 
wider  portion  of  the  test  section  width.  There  is  however,  a  reduction  in  the  amount  of 
mixed  fluid  (p)  in  the  middle  of  the  layer  as  the  entrainment  tongues  become  more 
prominent.  For  the  cases  forced  at/  =  2  Hz,  there  is  a  slight  variation  in  (p^)^  as  a 
function  of  x,  similar  to  the  unforced  layer.  However,  with  increasing  amplitude  there  is 
a  noted  decrease  in  (/?„,)  ,  see  Figures  14  -  16.  When  forced  at/=  4  Hz,  A  =  2°,  see 

Figure  17,  there  is  a  continual  decrease  in  the  amount  of  mixed  fluid  in  the  center  of  the 
layer.  At  the  larger  amplitudes  A  =  4°  and  6°  the  amount  of  mixed  fluid  continually 
decreases  except  at  x  =  20.0  cm  where  a  slight  increase  in  is  noted,  see  Figures 

18  and  19  (b).  The  cases  forced  at/=  8  Hz  show  markedly  different  behavior  than  the 
previous  cases.  When  forced  at  /  =  8  Hz,  A  =  2°  initially  decreases  then 

continually  increases;  while  (p„)^^  continually  increases  at  the  higher  amplitudes. 

The  reduction  in  the  amount  of  mixed  fluid  in  the  central  portion  of  the  layer  is 
accompanied  by  an  increase  in  the  amount  of  the  unmixed  fluid.  In  all  forced  cases  there 
is  more  pure  high-speed  fluid  than  pure  low-speed  fluid  in  the  center  of  the  layer.  This 
is  most  obvious  in  the  cases  forced  at/=  4,  8  Hz,  see  Figures  17  -  22  (c)  and  (d).  The 
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unequal  amounts  of  unmixed  fluid  in  the  center  of  the  layer  are  related  to  the  inherent 
asymmetry  of  the  shear  layer. 

It  was  previously  noted  in  the  flow  visualization  images  of  Figures  8-12  that 
some  of  the  forced  cases  were  influenced  by  the  finite  height  of  the  test  section.  This  is 
also  seen  in  the  total  mixed  fluid  probability  curves  in  Figures  16  (b)  (/■=  2  Hz,  A  =  8°) 

and  19  (b)  (/"=  4  Hz,  A  =  6°)  as  nonzero  values  of  at  y  =  -20  cm. 

\ 

3.2.2  Forcing  amplitude  dependence 

In  the  previous  section  the  streamwise  dependence  of  two-dimensional  forcing  was 
noted  for  the  various  forcing  conditions.  In  this  section  the  dependence  on  the  forcing 
amplitude  will  be  shown  at  a  fixed  downstream  location.  In  Figures  23  -  25,  transverse 
profiles  of  the  average  concentration  the  total  mixed  fluid  probability  (y)  and 
unmixed  fluid  probabilities,  pg  (y)  and  Pi  (y)  are  shown  as  a  function  of  the  forcing 
amplitude  atx  =  20.0  cm.  The  unforced  case  is  included  for  comparison.  Figures  for  the 
remaining  measurement  locations  are  provided  in  Appendix  A. 

The  'flat  region'  of  the  average  concentration  profiles  associated  with  the  passage 
of  well  mixed  vortex  cores  is  not  present  in  the  2  Hz  cases  because  the  vortices  are  not 
fully  formed  at  this  location.  There  is  almost  no  difference  in  the  average  concentration 
profiles  between  the  unforced  case  and  the  cases  forced  at  amplitudes  A  =  2°  and  4°,  see 
Figure  22  (a).  Well  mixed  vortex  cores  are  associated  with  the  well  defined  structures 
of  the  frequency  locked  region  seen  in  the  cases  forced  at  /  =  4  and  8  Hz.  The 
developing  'flat  region'  may  be  seen  in  the  4  Hz  cases  and  a  very  well  defined  'flat  region' 
is  observed  in  the/=  8  Hz  cases,  see  Figures  24  and  25  (a). 


(a)  Average  concentration 


(b)  Total  mixed  fluid  probability 


Figure  23.  Transverse  profiles  for  the  shear  layer  forced  at  /=  2  Hz  as  a  function  of 
amplitude  A,  at  v*  =  0.44  (x  =  20.0  cm),  compared  to  the  unforced  case. 
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(a)  Average  eoncentration 


(b)  Total  mixed  fluid  probability 


(e)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  24.  Transverse  profiles  for  the  shear  layer  forced  at/ =  4  Hz  as  a  function  of 
amplitude  A,  at  x*  =  0.88  (jc  =  20.0  cm),  compared  to  the  unforced  case. 
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(a)  Average  concentration 


y  (mm) 

(c)  Pure  low-speed  fluid  probability 
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(b)  Total  mixed  fluid  probability 
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(d)  Pure  high-speed  fluid  probability 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  25.  Transverse  profiles  for  the  shear  layer  forced  at/ =  8  Hz  as  a  function  of 
amplitude  A,  at  x*  =  1.78  (x  =  20.0  cm),  compared  to  the  unforced  case. 
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In  all  of  these  cases,  the  width  of  the  layer  is  larger  than  the  unforced  case  and 
increases  with  increasing  amplitude.  In  this  sense  the  effect  of  the  forcing  amplitude  on 
the  shear  layer  width  is  similar  to  increasing  the  downstream  distance.  The  cases  forced 
at  f  =  2  and  4  Hz  show  a  reduction  in  ( p  )  with  increasing  amplitude.  Coinciding 
with  the  reduction  in  the  amount  of  mixed  fluid  in  the  center  of  the  layer  there  is  an 
increase  in  the  amounts  of  pure  fluids  from  both  streams.  Note  that  there  is  in  general 
significantly  more  high-speed  fluid  in  the  center  and  on  the  lower  side  than  low-speed 
fluid  in  the  center  and  on  the  upper  side.  Again,  this  is  related  to  the  inherent 
entrainment  asymmetry  of  the  shear  layer.  The  cases  forced  at/=  8  Hz  exhibit  a  similar 
behavior.  Unlike  the  2  and  4  Hz  cases,  however,  there  is  no  dramatic  reduction  in 
(jPm)max  ^  function  of  amplitude  for  the  8  Hz  cases  at  this  x.  In  this  case,  is 

nearly  constant  at  this  location  (x  =  20.0  cm).  This  is  believed  to  be  connected  to  the 
saturation  of  the  large  structures  in  the  frequency  locked  region. 

3.3  Streamwise  variation  of  5.,  8^,  8^/8,  and  (p^) 

The  detailed  profiles  of  the  previous  section  can  be  used  to  compute  the  shear 
layer  width  and  other  mixing  quantities  such  as  the  mixed-fluid  thickness  8^ ,  mixed 
fluid  fraction  5  /S,  and  the  amount  of  fluid  in  the  center  of  the  layer  (p^)  .  These 

m  1  -I  \j.  m/jnsx 

results  are  discussed  in  this  section. 

In  Figures  26  (a)  -  (c)  the  shear  layer  width  Sj ,  is  plotted  versus  x  as  a  function 
of  amplitude  for  the  forcing  frequencies  /  =  2,  4  and  8  Hz,  respectively.  The  unforced 
layer  is  noted  to  have  linear  growth  rate.  From  this  Figure,  it  may  also  be  seen  that  the 
width  of  the  forced  layer  is  greater  than  the  unforced  layer  and  increases  with  downstream 
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X  (cm) 

o  Unforced  case,  o  A  =  2°,  v  A  =  4°,  •  A  =  6° 

Figure  26.  Streamwise  evolution  of  the  shear  layer  width  6i  as  a  function  of  the  forcing 
amplitude  A  compared  to  the  unforced  case. 


distance.  The  highest  amplitude  cases  all  seem  to  saturate  at  approximately  35  mm  in 
width.  The  forced  layer  is  expected  to  initially  have  enhanced  growth  rate,  then  taper  off 
(saturate)  and  finally  resume  the  natural  growth  rate,  see  section  1.2.  The  growth  of  the 
forced  layers  may  be  hindered  by  the  finite  size  (height)  of  the  test  section,  2h  =  40  mm, 

or  may  be  reaching  the  saturation  point  at  x*=  1. 

The  total  amount  of  mixed  fluid  in  terms  of  the  mixed  fluid  thickness  8^  is 
plotted  in  Figure  27.  The  amplitude  dependence  for  the  /  =  2  Hz  cases  are  shown  in 
Figure  27  (a).  Only  at  the  largest  amplitude  A  =  8°  was  there  a  significant  increase, 
approximately  11%,  in  the  mixed  fluid  thickness.  There  was  less  than  a  2%  change  in 


(mm) 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  27.  Streamwise  evolution  of  the  mixed  fluid  thickness  5„  as  a  function  of  the 
forcing  amplitude  A  compared  to  the  unforced  case. 
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o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 

Figure  27.  Streamwise  evolution  of  the  mixed  fluid  thickness  5„  as  a  function  of  the 
forcing  amplitude  A  compared  to  the  unforced  case. 


the  mixed  fluid  thickness  when  compared  to  the  unforced  shear  layer  at  the  two  smaller 
amplitudes.  Small  increases,  5%  at  most,  were  similarly  noted  for  the  cases  forced  at/ 
=  4  Hz,  see  Figure  27  (b).  There  was,  however,  a  marked  increase  in  the  amount  of 
mixed  fluid  at  all  amplitudes  when  the  layer  was  forced  at  8  Hz,  Figure  27  (c).  The 
largest  change  observed  in  this  study  was  29%,  when  the  layer  was  forced  at  /  =  8  Hz, 
A  =  6°. 

It  was  noted  above  that  the  largest  increase  in  the  mixed  fluid  thickness  compared 
to  the  unforced  layer  was  29%,  for  the  case  forced  at/  =  8  Hz,  A  =  6°.  The  shear  layer 
width  in  this  case  was  6i  =  35.2  mm,  an  increase  of  19%  over  the  unforced  layer  (5i  = 
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29.6  mm),  see  Figure  26  (c).  One  way  to  quantify  the  amount  of  mixing,  while  taking 
into  account  the  increasing  width  of  the  layer  is  to  use  the  mixed  fluid  fraction 
the  fraction  of  the  layer  width  occupied  by  mixed  fluid.  •  In  Figure  28,  the  mixed 

fluid  fraction  is  shown  as  a  function  of  amplitude.  In  most  cases  the  forced  layer  has  less 
mixed  fluid  per  unit  width  than  the  unforced  layer,  and  in  that  sense,  is  a  less  efficient 
mixer.  The  exceptions  are  the  cases  forced  at  /  =  8  Hz,  A  =  4°  and  6°  at  the  farthest 
downstream  locations  x  =  18.5  and  all  amplitudes  at  x  =  20.0  cm.  The  reduction  in  the 
mixed  fluid  fraction  as  a  result  of  forcing  was  similarly  reported  by  Koochesfahani  and 
McKinnon  (1991).  It  is  interesting  to  note  that  the  mixed  fluid  fraction  for  the  cases 
presented  here  are  at  best  no  greater  than  0.52;  only  slightly  more  than  half  of  the  fluid 


12.5  14.0  15.5  17.0  18.5  20.0 

X  (cm) 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 


Figure  28.  Streamwise  evolution  of  the  mixed  fluid  fraction  5„,/5i  as  a  function  of  the 
forcing  amplitude  A  compared  to  the  unforced  case. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 
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o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  28.  Streamwise  evolution  of  the  mixed  fluid  fraction  6„/6i  as  a  function  of  the 
forcing  amplitude  A  compared  to  the  unforced  case. 
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within  the  shear  layer  is  mixed  (to  within  the  measurement  resolution  200  ^m). 

The  amount  of  mixed  fluid  in  the  center  of  the  layer  is  characterized  by  the  peak 
of  the  mixed  fluid  probability  curve  The  amplitude  dependence  of  is 

shown  as  a  function  of  x  in  Figure  29.  It  is  observed  that  there  is  more  mixed  fluid  in 
the  center  of  the  unforced  shear  layer  compared  to  all  of  the  forced  cases.  For  the  cases 
forced  at/=  2  Hz,  (p^)  is  in  general  decreasing  with  amplitude  and  downstream 
distance.  The  4  Hz  cases  exhibit  a  similar  behavior  with  the  exception  of  a  slight 
increase  at  x  =  20.0  cm  for  the  two  larger  amplitudes.  When  forced  at/=  8  Hz,  A  =  2° 
(p  )  initially  decreases  then  continually  increases;  while  at  higher  amplitudes  the 
amount  of  mixed  fluid  in  the  center  of  the  layer  is  observed  to  continually  increase.  It 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 

Figure  29.  Streamwise  evolution  of  the  amount  of  mixed  fluid  in  the  center  of  the  layer 
as  a  function  of  the  forcing  amplitude  A  compared  to  the  unforced  case. 


12.5  14.0  15.5  17.0  18.5  20.0 

X  (cm) 

o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 

Figure  29.  Streamwise  evolution  of  the  amount  of  mixed  fluid  in  the  center  of  the  layer 
as  a  function  of  the  forcing  amplitude  A  compared  to  the  unforced  case. 
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is  interesting  to  note  that  reaches  a  minimum  (in  the  4  Hz  cases)  at  x  *  ~  1,  when 

the  structures  saturate;  this  will  be  further  discussed  in  section  3.4. 

3.4  Nondimensional  streamwise  variation  of  8^,  SJ8^  and 

The  previous  sections  discussed  the  evolution  of  the  shear  layer  in  terms  of  the 
downstream  coordinate  x;  now  flow  evolution  will  be  discussed  in  terms  of  the 
nondimensional  streamwise  coordinate  x*=  Xxf  I U^.  The  shear  layer  width  and 
mixed  fluid  thickness  8^  normalized  by  the  forcing  wavelength  UJf  are  shown  in 
Figure  30  and  31  versus  x*.  These  Figures  display  all  of  the  forced  cases  from  Figures 


* 

X 


o/=  2  Hz,  A  =  2°,  o/=  2  Hz,  A  =  4°,  v/=  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,A  =  6° 

+  /=  8  Hz,  A  =  2°,  ®/=  8  Hz,  A  =  4°,  s/=  8  Hz,  A  =  6° 

Figure  30.  Nondimensional  streamwise  evolution  of  the  shear  layer  width  normalized  by 
the  forcing  wavelength  for  the  various  forcing  conditions. 
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o/=  2  Hz,  A  =  2°,  □/=  2  Hz,  A  =  4°,  v/=:  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,  A  =  6° 

+  /=  8  Hz,  A  =  2°,  ®/=  8  Hz,  A  =  4°,  e/=  8  Hz,  A  =  6° 

Figure  3 1 .  Nondimensional  streamwise  evolution  of  the  mixed  fluid  thickness  normalized 
by  the  forcing  wavelength  for  the  various  forcing  conditions. 


26  and  27  in  nondimensional  form.  As  expected,  with  increasing  x*  the  shear  layer  width 
continually  increases.  The  saturation  of  the  structures  may  also  be  seen  in  the  frequency 
locked  region,  near  =  1.5.  The  amplitude  effects  are  also  observed  to  increase  the 
shear  layer  width.  Similarly,  the  amount  of  mixed  fluid  increases  with  increasing  x  * . 
There  is,  however,  no  saturation  effect  in  the  frequency  locked  region. 

The  mixed  fluid  fraction  8^/8^  is  displayed  in  Figure  32  as  a  function  of  x*. 
Noteworthy  features  of  this  Figure  are  a  local  minimum  occurring  near  x  *  =  0.8 ;  and  a 
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o/=  2  Hz,  A  =  2°,  □/=  2  Hz,  A  =  4°,  v/=  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,A  =  6° 

+  /=  8  Hz,  A  =  2°,  e /=  8  Hz,  A  =  4°,  s /=  8  Hz,  A  =  6° 

Figure  32.  Nondimensional  streamwise  evolution  of  the  mixed  fluid  fraction  for  the 

various  forcing  conditions. 


general  increase  in  the  mixed  fluid  fraction  for  x  *  >  1.0.  Figure  33  shows  the  probability 
of  mixed  fluid  in  the  center  of  the  layer  which  is  proportional  to  the  amount  of 

mixed  fluid  in  the  central  of  the  layer.  It  is  seen  that  ,  in  general,  decrease  while  x  * 

<1.0  and  increase  for  jc*  >  1.0.  These  observations  are  connected  to  the  growth  rate, 
size  and  spacing  of  the  structures.  The  local  minimum  appears  to  be  connected  to  the 
enhanced  growth  rate  of  the  layer;  in  this  region  the  layer  width  is  increasing  faster  than 
the  amount  of  mixed  fluid  in  the  layer.  In  the  frequency  locked  region  the  opposite  is 
true,  layer  growth  rate  is  not  increasing  as  fast  as  the  amount  of 
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o/=  2  Hz,  A  =  2°,  o/=  2  Hz,  A  =  4°,  v/=  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,A  =  6° 

+  /=8Hz,  A  =  2°,  ®/=8Hz,A  =  4°,  b/=8Hz,A  =  6° 

Figure  33.  Nondimensional  streamwise  evolution  of  the  amount  of  mixed  fluid  in  the 
center  of  the  layer  for  the  various  forcing  conditions. 


mixed  fluid  in  the  layer,  see  Figures  26  and  27.  Physically,  the  reason  that  the  frequency 
locked  region,  ie.  8  Hz  cases,  has  more  mixed  fluid  is  that  the  structures,  although 
roughly  the  same  size  as  the  4  Hz  cases,  are  twice  as  many  in  number. 

The  above  nondimensional  analysis  used  the  shear  layer  width  Sj  and  the  forcing 
wavelength  UJf  as  length  scales.  The  shear  layer  width,  recall  is  computed  similar  to 
boundary  layer  thickness,  in  that  the  edge  of  the  layer  is  defined  by  the  points  where  the 
total  mixed  fluid  probability  curve  reaches  1%  of  its  maximum  value.  It  suffers  from 
resolution  problems,  due  the  non-reacting  experiments  used  in  this  work  over-estimating 


59 


the  amount  of  mixed  fluid.  In  some  cases  the  p^(y)  curves  asymptotically  approach 
zero,  artificially  increasing  the  shear  layer  width.  Perhaps  a  more  robust  length  scale 
would  be  the  concentration  thickness  6^ ,  an  integral  quantity  analogous  to  momentum 
thickness  (introduced  in  section  2.4)  based  on  average  concentration  measurements.  The 
relationship  between  the  concentration  thickness  and  momentum  thickness  is  displayed 
in  Figure  34.  For  the  unforced  layer  the  relationship  is  nearly  linear,  as  the  frequency  and 
amplitude  increase  there  is  an  apparent  scatter  in  the  plot  at  large  .  This  scatter  may 
be  the  result  of  deficiencies  on  the  part  of  both  Sj  and  6^ .  The  concentration  thickness 
may  suffer  when  there  is  a  well  developed  flat  region  on  the  average  concentration 
profile. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 

Figure  34.  Streamwise  evolution  of  the  concentration  thickness  as  a  function  of  the 
forcing  amplitude  A  compared  to  the  unforced  case. 


o  Unforced  case,  □  A  =  2°,  ^  A  =  4°,  •  A  =  6° 
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o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  34.  Streamwise  evolution  of  the  concentration  thickness  as  a  function  of  the 
forcing  amplitude  A  compared  to  the  unforced  case. 
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The  concentration  thickness  normalized  by  the  forcing  wavelength  6^  //  is 
shown  versus  x  *  in  Figure  35;  note  the  similarity  to  Figure  30. 


o/=  2  Hz,  A  =  2°,  □/=  2  Hz,  A  =  4°,  v/=  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,A  =  6° 

+  /=  8  Hz,  A  =  2°,  ®/=  8  Hz,  A  =  4°,  ®/=  8  Hz,  A  =  6° 

Figure  35.  Nondimensional  streamwise  evolution  of  the  concentration  thickness 
normalized  by  the  forcing  wavelength  for  the  various  forcing  conditions. 


The  mixed  fluid  thickness  normalized  by  the  concentration  thickness  8^16^,  Figure  36, 
shows  a  trend  very  similar  to  the  mixed  fluid  fraction  versus  x*,  see  Figure  32.  The 
most  notable  feature,  again,  is  the  local  minimum  occurring  in  the  vicinity  of  x*  =  1.0, 
the  end  of  the  enhanced  growth  rate  region.  As  x*  increases  past  1.0,  a  steady  increase 
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in  8^  16^  is  observed.  Trends  in  Figure  36  are  the  same  as  those  shown  in  Figure  32,  the 
only  difference  being  the  lack  of  the  scatter  in  Figure  32. 


o/=  2  Hz,  A  =  2°,  □/=  2Hz,  A  =  4°,  v/=  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,A  =  6° 

+  /=  8  Hz,  A  =  2°,  ®/=  8  Hz,  A  =  4°,  s/=  8  Hz,  A  =  6° 

Figure  36.  Nondimensional  streamwise  evolution  of  the  mixed  fluid  thickness  normalized 
by  the  concentration  thickness  for  the  various  forcing  conditions. 


The  plots  in  this  section  were  normalized  with  respect  to  the  forcing  frequency 
only;  effects  of  the  forcing  amplitude  are  not  effected  by  this  normalization.  It  is  not 
clear  at  this  time  if  the  amplitude  effect  can  be  properly  normalized. 


3.5  Probability  density  function  of  ^ 
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Results  thus  far  have  addressed  mixing  in  terms  of  mixed  fluid  at  all 
concentrations.  The  influence  of  forcing  on  the  actual  distribution  of  the  mixed  fluid 
composition  will  now  be  discussed  in  terms  of  the  average  probability  density  function 
P  (^,  y).  According  to  Koochesfahani  and  MacKinnon  (1991),  the  use  of  a  single  pdf  to 
characterize  the  composition  field  of  both  the  unforced  and  the  forced  layers  is  warranted 
here,  since  the  shape  of  the  pdf  was  observed  to  be  nearly  invariant  across  the  layer 
width.  The  composition  distribution  for  the  unforced  and  forced  shear  layers  are 
presented  in  Appendix  B. 

3.5.1  Streamwise  variation  of  the  total  pdf 

Figures  37  displays  the  streamwise  evolution  of  the  total  pdf  of  ^  for  the  unforced 
shear  layer.  The  general  shape  of  the  pdf  is  preserved  as  downstream  distance  increases. 
An  increase  in  mixed  fluid  with  increasing  downstream  distance  is  seen  by  the  increased 
height  of  the  peak  of  the  pdf  located  near  ^  =  0.75.  This  peak  in  the  pdf  will  be  referred 
to  as  the  most  probable  concentration.  It  is  interesting  to  note  how  the  pdf  fills  in  around 
the  most  probable  concentration,  but  not  at  very  low  concentrations.  It  is  also  observed 
that  the  peak  tends  to  shift  slightly  toward  lower  values  as  x  increases.  Initially  the  layer 
has  an  excess  of  high-speed  fluid  (^  =  1)  resulting  in  a  most  probable  concentration  at 
very  high  values  of  As  the  mixing  process  continues  the  most  probable  concentration 
shifts  toward  the  low-speed  side  =  0).  This  result  was  first  reported  by  Koochesfahani 
and  Dimotakis  (1986). 

The  streamwise  dependence  of  the  pdf  for  the  shear  layer  forced  at  /  =  2,  4  and 
8  Hz,  A  =  4°  are  displayed  in  Figures  38  -  40,  respectively.  The  remaining  forced  cases 
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are  presented  in  Appendix  C.  In  all  cases  mixing  is  enhanced  mostly  at  high  values  of 
concentration  with  increasing  downstream  distance.  Unlike  the  unforced  layer  which 
shows  a  slight  shift  in  the  most  probable  concentration  to  lower  values,  some  of  the  more 
highly  forced  cases  exhibit  a  slight  shift  toward  higher  values,  see  Figures  39  and  40. 


o  X  =  12.5,  □  jc  =  14.0,  ^  x  =  15.5,  •  x  =  17.0,  ■  x  =  18.5,  a  x  =  20.0  cm 
Figure  37.  Streamwise  evolution  of  the  total  pdf  of  the  unforced  shear  layer. 


o  JC  =  0.28,  □  X  =  0.31,  V  X  =  0.34,  •  x  =  0.37,  ■  x  =  0.41,  a  x  =  0.44 

Figure  38.  Nondimensional  streamwise  evolution  of  the  total  pdf  for  the  shear  layer 
forced  at  /  =  2  Hz,  A  =  4°. 


Figure  39.  Nondimensional  streamwise  evolution  of  the  total  pdf  for  the  shear  layer 
forced  at/  =  4  Hz,  A  =  4°. 
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o  /  =  1.11,  □  /  =  1.24,  V  X*  =  1.38,  •  X*  =  1.51,  m  X*  =  1.64,  a  /  =  1.78 

Figure  40.  Nondimensional  streamwise  evolution  of  the  total  pdf  for  the  shear  layer 
forced  at  /  =  8  Hz,  A  =  4°. 

3.5.2  Amplitude  dependence  of  the  total  pdf 

The  forcing  amplitude  dependence  is  shown  in  Figures  41-43  for  the  cases  forced 
at  /  =  2,  4  and  8  Hz  at  the  farthest  downstream  location  x  =  20.0  cm,  corresponding  to 
X*  =  0.44,  0.88,  1.78,  respectively.  The  remaining  cases  are  presented  in  Appendix  D. 
The  main  result  is  that  in  all  cases  mixing  is  enhanced  mostly  at  high  values  of 
concentration.  The  effect  of  the  forcing  amplitude  is  simply  to  increase  the  proportion 
of  fluid  mixed  at  higher  values  of  concentration.  In  addition  to  enhanced  mixing,  there 
is  a  change  in  the  composition  field.  In  all  cases  the  most  probable  concentration  is 
shifted  slightly  toward  higher  values  of  A  shift  in  the  predominate  mixed-fluid 
concentration  was  also  noted  by  Koochesfahani  and  M“‘'Kinnon  (1991). 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8 


Figure  41.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/=  2  Hz, 
at  X  =  20.0  cm  (x*  =  0.44),  compared  to  the  unforced  case. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 

Figure  42.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/ =  4  Hz 
at  X  =  20.0  cm  (x*  =  0.88),  compared  to  the  unforced  case. 
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o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  43.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/=  8  Hz, 
at  x  =  20.0  cm  (x*  =  1.78),  compared  to  the  unforced  case. 


3.6  Comparison  with  Koochesfahani  and  M^'Kinnon  (1991) 

The  current  flow  conditions,  U,  =  40  cm/sec  and  U2  =  20  cm/sec,  were  chosen  to 
be  the  same  as  those  in  Koochesfahani  and  M‘“^Kinnon  (1991),  to  allow  comparison  of  the 
forcing  mechanism  on  the  development  of  the  forced  shear  layer.  The  facility  used  in  the 
current  experiment  was  the  same  facility  used  by  Koochesfahani  and  M‘“'Kinnon  (1991). 
In  the  previous  work,  data  were  presented  at  a  single  downstream  position,  x  =  20  cm, 
where  the  local  Reynolds  number  Reg  =  Af/Sj/v,  was  approximately  6,600.  In  the 
current  work  Re.  =  6,000  at  x  =  20.0  cm.  Due  to  the  finite  size  (height  and  length)  of 
the  test  section  it  was  not  possible  to  match  the  Reynolds  numbers  of  the  two 
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experiments. 

As  stated  earlier,  the  unforced  layer  in  the  current  study  is  not  the  same  as  the 
natural  layer,  due  to  the  presence  of  the  airfoil;  the  main  difference  being  a  reduction  in 
the  width  of  the  unforced  layer.  A  summary  of  shear  layer  width  5i  ,  mixed  fluid 
thickness  and  the  mixed  fluid  fraction  SJS^,  for  the  natural  and  unforced  layers  are 
given  in  Table  2. 

Table  2.  Summary  of  shear  layer  width,  mixed  fluid  thicknesses  and 
mixed  fluid  fraction  for  the  natural  and  unforced  layers. 


(mm) 

5i  (mm) 

Natural 

16.1 

33.3 

0.483 

Unforced 

14.4 

29.6 

0.486 

The  natural  layer  of  Koochesfahani  and  ]VT'''Kinnon  has  width  8,  =  33.3  mm  and  the 
unforced  layer  in  the  current  investigation  was  bj  =  29.6  cm.  This  is  approximately  an 
11%  reduction.  The  reduction  in  the  width  of  the  layer  due  to  the  presence  of  an  airfoil 
is  consistent  with  the  results  reported  by  Koochesfahani  and  Dimotakis  (1989).  The  peak 
of  the  mixed  fluid  probability  curve,  in  the  central  region  of  the  shear  layer  is 

nearly  the  same,  implying  that  the  amount  of  mixed  fluid  atx  =  20.0  cm  is  not  changed 
due  to  the  presence  of  the  airfoil.  Therefore,  the  changes  in  between  the  unforced  and 
natural  layers  are  due  purely  to  the  reduced  width  of  the  shear  layer. 

In  the  forced  layers  of  Koochesfahani  and  M‘*'=Kinnon  (1991),  perturbations  were 
induced  by  sinusoidally  oscillating  the  high-speed  freestream  at  frequencies  of  4  and  8 
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Hz,  with  rms  amplitude  of  about  1.0%  of  the  freestream  speed.  There  is  some  ambiguity 
in  the  comparing  these  cases  to  the  present  work  since  the  perturbation  methods  are 
different  and  the  amplitudes  are  characterized  differently.  It  is  not  clear  if  the  methods 
by  which  each  of  these  forcing  techniques  are  quantified  are  related  to  each  other.  The 
effect  of  the  different  forcing  mechanisms  on  the  scalar  mixing  field  does,  however,  show 
the  same  trend.  In  Figure  44,  the  forcing  frequency  and  amplitude  dependence  of  the 


f(Hz) 


* 

X 


Figure  44.  Nondimensional  streamwise  (frequency)  dependence  of  the  mixed  fluid 
fraction  for  various  forced  cases  compared  with  Koochesfahani  and 
M^'^Kinnon  (1991). 
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mixed-fluid  fraction  for  the  two  forcing  mechanisms  are  displayed.  At  a  fixed 
downstream  location  the  nondimensional  streamwise  coordinate  x*  =  Xxf  I is  directly 
proportional  to  the  forcing  frequency  /.  Therefore  the  frequency  scale  of  Figure  44  may 
also  be  viewed  as  an  evolutionary  scale.  Only  for  the  cases  forced  at/ =  8  Hz  is  there 
more  mixed-fluid  per  unit  width  of  the  layer  than  in  the  unforced  case  at  this  downstream 
location;  if  the  flow  is  examined  farther  downstream,  further  increases  in  the  mixing 
would  be  expected.  Table  3  summarizes  the  shear  layer  width  8i ,  mixed  fluid  thickness 
6^  and  the  nnixed  fluid  fraction  SJS^  for  the  forced  cases  in  both  of  these  studies. 

Table  3.  Summary  of  shear  layer  width,  mixed  fluid  thickness  and  mixed 
fluid  fraction  for  the  different  forcing  conditions. 


(mm) 

6i  (mm) 

Koochesfahani  and  M“‘^Kinnon 

/=4Hz 

18.3 

40.5 

0.452 

/=8Hz 

19.3 

37.5 

0.515 

Present  work 

/=  2  Hz,  A  =  2° 

14.2 

29.4 

0.482 

/=2Hz,  A  =  4° 

14.4 

31.0 

0.463 

/=  2  Hz,  A  =  8° 

15.7 

34.7 

0.453 

/=4Hz,  A  =  2° 

14.5 

30.6 

0.472 

/=  4  Hz,  A  =  4° 

14.8 

34.2 

0.433 

/=4Hz,  A  =  6° 

15.6 

36.0 

0.434 

/=  8  Hz,  A  =  2° 

15.5 

31.2 

0.498 

/=8Hz,  A  =  4° 

17.6 

34.0 

0.517 

/=  8  Hz,  A  =  6° 

18.5 

35.2 

0.526 
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The  total  pdf  for  the  oscillating  freestream  experiments  showed  trends  similar  to 
the  present  work.  In  both  cases  the  effect  of  forcing  was  to  increase  mixing  at  mainly 
higher  values  of  concentration.  It  is  noted,  however,  that  the  predominant  concentration 
of  mixed  fluid  shows  only  slight  variations  relative  to  the  unforced  case,  in  contrast  to 
the  oscillating  freestream  experiment  which  showed  a  larger  shift  to  higher  concentrations. 
It  is  not  clear  at  this  time  why  the  current  forcing  method  results  in  smaller  shifts  of  the 
predominant  mixed  fluid  concentration.  It  is  possible  that  the  wake  of  the  airfoil  may  be 
a  contributing  factor. 


3.7  Comparison  with  Roberts  and  Roshko  (1985) 

In  the  work  of  Roberts  and  Roshko  (1985)  the  effect  of  a  periodic  disturbance  on 
the  amount  of  chemical  product  in  a  forced  shear  layer  was  examined  at  pre-transitional 
and  post-transitional  Reynolds  numbers.  Forcing  was  achieved  by  sinusoidally  oscillating 
the  high-speed  freestream,  at  2%  (rms)  of  the  freestream  speed.  The  chemical  product 
was  measured  using  a  laser  absorption  technique.  In  this  work  the  amount  of  mixing  was 
quantified  by  the  product  thickness,  defined  by 


where  C-  is  the  low-speed  feestream  reactant  concentration  and  c  (y)  is  the  average 

Zq  y 

product  concentration  defined  as 


1-e 


e 
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The  product  thickness  for  the  various  forcing  conditions  was  noted  relative  to  the  product 
thickness  of  the  unforced  layer,  denoted  by  the  ratio  where  8^^  is  the 

product  thickness  and  the  subscripts  F  and  N  refer  to  the  forced  and  natural  layers, 
respectively. 

In  these  experiments  the  unit  Reynolds  number  Re '  =  At//  v  of  the  pre-transitional 
(low  Reynolds  number)  case  was  905  cm  Vsec,  and  for  the  post-transitional  (high 
Reynolds  number)  case  4340  cm  Vsec.  The  natural  frequency  for  the  low  Reynolds 
number  case  was/g  =  8.3  Hz  and  for  the  high  Reynolds  number  case/g  =  62.7  Hz.  In  the 
current  work  the  unit  Reynolds  number  is  2000  cm  Vsec,  corresponding  to  the  later  stages 
of  the  mixing  transition.  The  natural  frequency  for  the  unforced  case  is  /g  ~  27  Hz. 

It  was  found,  in  the  work  of  Roberts  and  Roshko,  for  the  low  Reynolds  number 
case  that  forcing  had  steadily  increased  the  ratio  of  the  product  thicknesses  {pp)Fli^p)N 
through  the  enhanced  growth  region,  reaching  a  local  maximum  near  x*=2.  Once  past 
the  frequency  locked  region  the  product  thicknesses  steadily  decreased.  The  ratio  of 
forced  product  thickness  to  the  natural  product  thickness  {8 is  shown  in  Figure 

45  as  a  function  of  x  *  for  the  low  Reynolds  number  case.  In  the  high  Reynolds  number 
case,  the  forced  product  thickness  is  initially  as  much  as  4  times  larger  than  the  unforced 
product  thickness;  their  ratio  (5^^)^ steadily  decreases  through  the  enhanced  growth 
region.  In  the  frequency  locked  region  and  beyond,  x*>l,  the  ratio  of  the  forced  to  the 
natural  product  thickness  remains  slightly  less  than  the  one,  see  Figure  46.  This  indicates 
that  there  is  actually  slightly  less  chemical  product,  on  average,  in  the  forced  layer  than 
in  the  unforced  layer. 

In  the  present  work,  the  ratio  of  the  forced  mixed  fluid  thickness  to  the  unforced 
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Figure  45.  Ratio  of  forced  product  thickness  to  natural  product  thickness  as  a  function 
of  X*  :  low  Reynolds  number  case,  reprinted  from  Roberts  (1985). 


Figure  46.  Ratio  of  forced  product  thickness  to  natural  product  thickness  as  a  function 
of  /  :  high  Reynolds  number  case,  reprinted  from  Roberts  (1985). 
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mixed  fluid  thickness  is  analogous  to  the  ratio  of  the  forced  to  natural 

product  thickness,  see  Figure  47.  It  should  be  noted  that  the  current  experiments  were 
done  using  a  non-reacting  technique  which  over-estimates  the  amount  of  mixed  fluid. 
Therefore,  the  trends  shown  in  the  data  are  of  interest  not  the  actual  values  of  the 
quantities.  The  trends  observed  in  the  present  data  are  most  similar  to  the  trends  shown 
in  the  low  Reynolds  number  case.  Figure  45.  In  the  enhanced  growth  region  there  is  a 
slight  increase  in  mixed  fluid  thickness  over  the  unforced  case;  and  in  the  frequency 
locked  region  there  is  a  large  increase  in  8^  that  seems  to  peak  near  x*=  2.  Note  the 
expanded  scale  of  Figure  47. 


o/=  2  Hz,  A  =  2°,  □/=  2  Hz,  A  =  4°,  v/=  2  Hz,  A  =  8° 

•  /=4Hz,A  =  2°,  ■/=4Hz,A  =  4°,  a/=4Hz,A  =  6° 

+  /=  8  Hz,  A  =  2°,  ®/=  8  Hz,  A  =  4°,  ^/=  8  Hz,  A  =  6° 

Figure  47.  Ratio  of  forced  product  thickness  to  unforced  product  thickness  as  a  function 

r% 

of  X  . 


Chapter  4 

CONCLUSIONS 


The  structure  and  the  scalar  mixing  field  in  a  shear  layer  forced  by  an  oscillating 
airfoil  were  investigated  over  a  range  of  amplitudes,  frequencies  and  streamwise  locations 
using  laser  induced  fluorescence  in  the  passive  scalar  mode.  One  purpose  of  this  study 
was  to  compare  the  mixing  characteristics  of  a  shear  layer  which  is  forced  by  two 
different  methods.  Comparisons  were  made  with  previous  results  from  a  shear  layer 
perturbed  by  the  oscillation  of  one  freestream. 

It  was  observed  that  the  general  structure  and  the  growth  characteristics  of  the 
forced  shear  layer  here  were  similar  to  previous  results.  As  either  the  forcing  frequency 
or  amplitude  was  increased,  there  was,  in  most  cases,  a  significant  increase  in  the  shear 
layer  width.  Results  show  a  reduction  of  the  amount  of  mixed  fluid  in  the  central  portion 
of  the  forced  layer,  an  increase  in  the  width  of  the  region  where  mixed  fluid  is  found,  and 
an  increase  in  the  total  amount  of  mixed  fluid  in  the  layer.  The  fraction  of  the  layer 
width  occupied  by  mixed  fluid  was  found  to  vary  slightly,  however,  when  compared  to 
the  unforced  case.  These  results  are  consistent  with  those  from  a  shear  layer  forced  by 
the  oscillating  freestream  method.  The  predonoinant  mixed  fluid  concentration  showed 
similar  variation  as  a  result  of  forcing,  when  compared  to  previous  results  from  oscillating 
freestream  forcing  studies.  The  composition  distribution  for  the  unforced  and  the  forced 
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shear  layers  was  reported  to  be  essentially  uniform  across  the  width  of  the  layer  again. 
This  result  is  consistent  with  the  oscillating  freestream  experiment. 

The  effect  of  downstream  evolution  was  found  to  increase  shear  layer  width  and 
increase  the  total  amount  of  mixed  fluid  in  the  layer.  The  amount  of  rmxed  fluid  in  the 
center  of  the  layer  and  the  mixed  fluid  thickness  were  found  to  vary  depending  the  stage 
of  the  evolution  as  defined  by  the  nondimensional  streamwise  parameter  x  -  Xxf  /U^. 
A  local  minimum  was  reported  near  the  end  of  the  enhanced  growth  region  for  both  the 
mixed  fluid  fraction  and  the  amount  of  mixed  fluid  in  the  center  of  the  layer.  The  choice 
of  length  scale  when  computing  the  mixed  fluid  thickness  as  noted  to  have  a  slight  effect 
on  the  reported  efficiency  of  the  mixing  layer.  In  most  cases  the  unforced  layer  was 
observed  to  be  more  efficient,  in  terms  of  the  amount  of  mixed  fluid  per  unit  width  of  the 
layer,  than  the  forced  layer. 


APPENDICES 


APPENDIX  A 


Amplitude  dependence  of  the  transverse  profiles 


Amplitude  dependence  of  the  transverse  profiles  of  average  concentration  total 
mixed  fluid  probability  p^,  pure  low-speed  fluid  probability  and  pure  high-speed  fluid 
probability  p^  at  the  downstream  locations  not  presented  in  Section  3.2.2  are  shown  in 
Figures  A.l  -  A.  15. 
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(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


Figure  A.l. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 


Transverse  profiles  for  the  shear  layer  forced  at/ =  2  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  x*  =  0.28,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


Figure  A.2. 


o  Unforced  case,  □  A  =  2°,  ^  A  =  4°,  •  A  =  8° 


Transverse  profiles  for  the  shear  layer  forced  at/ =  2  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  jc*  =  0.31,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


Figure  A.  3. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 


Transverse  profiles  for  the  shear  layer  forced  at  /=  2  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  /  =  0.34,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 


Figure  A.4.  Transverse  profiles  for  the  shear  layer  forced  at  /=  2  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  x  =  0.37,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


Figure  A.5. 


o  Unforced  case,  o  A  =  2°,  v  A  =  4°,  •  A  =  8° 


Transverse  profiles  for  the  shear  layer  forced  at/ =  2  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  /  =  0.41,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


Figure  A.6. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Transverse  profiles  for  the  shear  layer  forced  at/  =  4  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  x  =  0.56,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


Figure  A.7. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Transverse  profiles  for  the  shear  layer  forced  at  /=  4  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  jc*  =  0.62,  compared  to  the  unforced  case. 


(a)  Average  concentration 


(b)  Total  mixed  fluid  probability 
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(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  A.ll.  Transverse  profiles  for  the  shear  layer  forced  at/ =  8  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  jc*  =  1.11,  compared  to  the  unforced  case. 


(c)  Pure  low-speed  fluid  probability 


(d)  Pure  high-speed  fluid  probability 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  A.  12.  Transverse  profiles  for  the  shear  layer  forced  at  /=  8  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  x*  =  1.24,  compared  to  the  unforced  case. 


(a)  Average  concentration 


(b)  Total  mixed  fluid  probability 


(c)  Pure  low-speed  fluid  probability  (d)  Pure  high-speed  fluid  probability 

o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  A.  15.  Transverse  profiles  for  the  shear  layer  forced  at/  =  8  Hz  as  a  function  of 
the  forcing  amplitude  A,  at  /  =  1.64,  compared  to  the  unforced  case. 


APPENDIX  B 


Distribution  of  mixed  fluid  composition 


The  distribution  of  mixed  fluid  composition  across  the  height  of  the  test  section 
for  the  unforced  and  forced  shear  layers  are  shown  in  Figures  B.l  -  B.IO  as  a  function 
of  the  streamwise  coordinate.  From  these  Figures  it  is  noted  that  the  pdf  is  essentially 
uniform  across  the  height  of  the  test  section,  similar  to  the  results  of  Koochesfahani 
(1984)  and  Koochesfahani  and  M^‘=Kinnon  (1991). 
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APPENDIX  C 


Streamwise  evolution  of  the  total  pdf 


Streamwise  evolution  of  the  total  pdf  for  the  forced  shear  layers  not  presented  in 
Section  3.5.1  are  shown  in  Figures  C.l  -  C.6. 
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APPENDIX  D 


Amplitude  dependence  of  the  total  pdf 

Amplitude  dependence  of  the  total  pdf  for  the  forced  shear  layers  at  the 
downstream  locations  not  presented  in  Section  3.5.2  are  shown  in  Figures  D.l  -  D.15. 
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o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 

Figure  D.l.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/  =  2 
Hz,  at  ;c  =  12.5  cm  (x*  =  0.28),  compared  to  the  unforced  case. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8° 

Figure  D.2.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/=  2 
Hz,  at  X  =  14.0  cm  (jc*  =  0.31),  compared  to  the  unforced  case. 


Figure  D.4.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at  /  =  2 
Hz,  at  X  =  17.0  cm  (x*  =  0.37),  compared  to  the  unforced  case. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  8 


Figure  D.5.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  foreed  at/  =  2 
Hz,  at  x=  18.5  em  (x*  =  0.41),  compared  to  the  unforced  ease. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 


Figure  D.6.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/  =  4 
Hz,  at  X  =  12.5  cm  (x*  =  0.56),  compared  to  the  unforced  case. 


Figure  D.IO.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/=  4 
Hz,  at  X  =  18.5  cm  (x*  =  0.82),  compared  to  the  unforced  case. 


Figure  D.12.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/=  8 
Hz,  at  jc  =  14.0  cm  (;c*  =  1.24),  compared  to  the  unforced  case. 


o  Unforced  case,  □  A  =  2°,  v  A  =  4°,  •  A  =  6° 

Figure  D.14.  Amplitude  dependence  of  the  total  pdf  for  the  shear  layer  forced  at/=  8 
Hz,  at  X  =  17.0  cm  (x*  =  1.51),  compared  to  the  unforced  case. 
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